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Abstract 


Great  effort  is  being  focused  on  making  the  next  generation  of  naval  combatant  ships 
more  resistant  to  the  effects  of  close-aboard  explosions.  The  examination  of  the 
deformation  modes  in  blast-loaded  metal  plating  suggests  that  a  physical  model  can  be 
developed  to  simulate  the  force  vs.  displacement  history  produced  by  an  impinging  shock 
wave  during  the  holing  phase.  Similar  approaches  have  been  successfully  used  to 
approximate  damage  due  to  grounding  and  ballistic  penetrators. 

In  this  case,  the  deformation  of  the  clamped  plate  is  modeled  in  two  stages:  (1)  dishing, 
which  leads  to  disking  and  (2)  radial  crack  propagation,  which  results  in  petalling.  In  the  first 
stage,  a  thin  geometrically-scaled  (0.90  mm,  1.15  mm,  and  1.40  mm  thick  by  300  mm  square) 
mild  steel  sheets  are  dished  inward  using  spherical  indenters  of  radii  20  mm,  50mm,  and  75 
mm.  The  sheets  have  an  average  tensile  strength  of  317  MPa  and  a  Rockwell  Superficial 
Hardness  Number  of  72  (HN15T72). 

This  portion  of  the  test  approximates  the  initial  material  stretching  done  by  a  spherical 
wave  at  various  standoff  distances.  The  spherical  indenter  produces  a  circular  hole,  which 
simulates  the  disk  of  material  normally  ejected  as  a  blast  front  penetrates  a  plate  section.  As 
the  material  reaches  a  critical  necking  thickness  at  the  edges  of  the  hole,  radial  cracks  form 
creating  petals.  During  the  second  stage,  an  oblique  conical  punch  is  used  to  simulate  the 
expanding  wave  front,  which  drives  open  the  petals,  causing  the  cracks  to  propagate  towards 
the  plate’s  clamped  boundaries.  By  measuring  the  resultant  forces  and  minimizing  the 
effects  of  friction,  the  total  bending  and  membrane  work  can  be  reasonably  estimated. 
Ultimately,  the  approximate  blast  damage  for  a  given  ship’s  hull  may  be  related  to  a  given 
charge  size  and  standoff  distance. 
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Introduction 


In  response  to  the  growing  number  of  threat  nations  and  the  increased  proliferation  of 
anti-ship  weaponry,  the  Office  of  Naval  Research  (ONR)  has  prompted  new  research  in  the 
areas  of  weapons  effects  and  ship  vulnerability.  Mindful  of  budgetary  pressures,  the  United 
States  Navy  is  exploring  the  most  cost-effective  methods  of  increasing  the  damage  resistance 
and  improving  the  overall  battle  effectiveness  of  its  warships.  (Refer  to  Appendix  A  for 
more  detailed  prefatory  material.) 

To  date,  the  plastic  deformation  resulting  from  a  close-proximity  explosion  (either  above 
or  below  the  waterline)  has  not  been  the  focus  of  any  large-scale,  publicly-accessible  research 
project.  Many  organizations  have  intensely  researched  topics  which  apply  similar  mechanics, 
such  as  ballistic  penetration,  collision  and  grounding  damage,  and  axial  tube  splitting.  There 
has  been  preliminary  work  by  both  Wierzbicki  (1996)/ (1999)  and  Nurick  (1996),  which 
examined  plate  tearing  and  petalling  done  by  on  contact  explosives.  (A  detailed  review  of 
past  work  in  this  area  and  the  associated  literature  can  be  found  in  Appendix  B) 

In  keeping  with  the  goals  of  the  U.S.  Navy  and  ONR,  this  research  focuses  on  the 
deformation  and  fracture  of  hull  plating,  which  is  subjected  to  either  an  underwater 
(UNDEX)  or  air  explosion.  Ultimately,  the  objective  of  the  study  is  to  provide  a  simplified 
method  of  gathering  benchmark  data  to  quantify  a  material’s  sensitivity  to  explosive  damage. 
An  inexpensive  two-stage  quasi-static  indentation  test  is  used  to  model  the  force  vs. 
displacement  history  of  a  thin  clamped  steel  sheet.  These  experimental  results  are  then 
compared  to  both  approximate  analytical  solutions  and  results  obtained  through  Numerical 
simulations  created  in  ABAQUS  (Static  Load  with  Pressure  Boundary  Conditions)  and  LS- 
DYNA  (Quasi-static  Load  with  Displacement  Boundary  Conditions). 
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Nomenclature 


R  Plate  Width 

r  Hole  Radius 

Rt,  Spherical  Punch  Radius 

Rc  Clamped  Plate  Width 

Ff  Frictional  Force 

P  Indenter  Force 

V  Indenter  Velocity 

p  Frictional  Coefficient 

oy  Material  Yield  Strength 

%  Material  Ultimate  Strength 

aQ  Material  Flow  Stress 

8  Plate  Central  Displacement 

6  Cone  Punch  Angle 

<pc  Indenter  Wrapping  Angle 

t  Plate  Thickness 

E  Indenter  Work 

1  Petal  Length 


Statement  of  Problem 


To  improve  blast  damage  resistance  in  future  classes  of  combatant  ships,  the  U.  S.  Navy 
and  ONR  are  entertaining  new  concepts  in  both  design  and  materials.  Although,  the 
characteristics  of  traditional  mild  and  high  strength  steels  have  been  studied  extensively  in 
naval  applications,  no  simple,  reliable  method  exists  to  predict  hull  panel  blast  damage  due  to 
a  close  proximity  blast  load.  Blast  holing  prediction  is  a  key  factor  in  assessing  ship 
survivability,  including  the  number  of  flooded  compartments,  the  ships  residual  section 
modulus,  and  probability  of  recovery.  Blast  damage  computer  codes  are  plentiful,  but  the 
accuracy  of  their  results  often  remains  in  question.  Consequently,  nearly  all  trusted  damage 
prediction  is  done  through  scaled  live  fire  testing,  which  is  both  expensive  and  time 
consuming. 


DISHED  SECTION 


Figure  1:  Simplified  Blast  Loaded  Plate  Geometry 
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The  long-term  challenge  of  this  research  is  to  combine  theories  from  tearing  fracture, 
plate  cutting,  petalling,  and  blast  loading  so  as  to  provide  a  preliminary  estimate  of  the 
damage  caused  by  a  given  charge  at  some  small  distance  from  a  clamped  metal  plate. 

The  simplified  physical  model  in  Figure  1,  presents  the  geometiy  of  the  problem.  Figure  2 
shows  clear  deformation  and  fracture  similarities  between  dynamic  (Rajendran  (2000))  and 
quasi-static  test  specimens. 


On  Contact  Blast  Test  (20  g  PEK-1  explosive)  Quasi-statically  Indented  Plate 


Figure  2:  Dynamically  and  Quasi-statically  Fractured  Specimens 

The  quasi-static  tests  have  two  immediate  purposes.  First,  if  the  time-pressure  history 
of  an  explosive  can  be  related  to  the  force  displacement  history  .of  the  proposed  quasi-static 
tests,  it  becomes  possible  to  relate  the  energy  of  the  incident  blast  wave  to  the  bending, 
stretching,  and  tearing  work  done  in  the  material.  Secondly,  by  examining  crack  initiation 
and  propagation  and  measuring  incremental  strains  during  these  quasi-static  tests,  one  can 
develop  a  rudimentary  crack-propagation  criterion.  Such  a  criterion  could  then  be  used  to 
improve  existing  finite  element  codes  and  coupled  hydro-codes.  It  is  believed  that  these 
tools  will  lead  to  improve  damage  estimates  (e.g.  resultant  hole  size  for  a  given  charge).  If 


such  estimates  can  be  made  quickly  and  inexpensively,  research  facilities  could  rapidly 
examine  the  suitability  of  a  wide  range  of  materials  and  structural  arrangements  with  minimal 


Formulation  of  Problem 

Experimental  Approach 

A  series  of  penetrator  tests  is  used  to  model  the  deformation  which  occurs  in  mild  steel 
panels  due  to  close  proximity  explosions.  The  modeling  occurs  in  two  stages.  During  the 
first  stage,  forces,  displacements,  and  local  strains  are  measured  as  a  hemispherical  punch  is 
used  to  dish  the  center  of  a  clamped  thin  mild  steel  sheet.  As  the  material  stretches  and 
thins,  a  critical  necking  thickness  will  be  reached  near  the  center  of  the  sheet.  A 
circumferential  crack  then  forms  at  the  necked  location.  Radial  cracks  will  subsequently  form 

at  the  newly  formed  hole’s  edge.  The  initiation  of  radial  cracks  will  mark  the  conclusion  of 
the  first  stage. 

In  the  second  stage,  an  oblique  conical  punch  is  used  to  examine  the  behavior  of  the 
material  while  petalling.  The  cone  is  used  to  propagate  either  pre-cut  or  naturally  formed 
radial  cracks,  causing  several  distinct  outward-opening  petals  to  form.  With  Teflon-based 
dry  film  lubricants  baked  onto  both  the  test  specimens  and  indenter  contact  surfaces, 
friction  is  kept  to  a  minimum  during  the  test.  Fine  grid  markings  (2.5  mm  squares)  are 
made  on  the  plate  in  order  to  compute  final  local  strains  at  periodic  locations  along  the 
plate  s  radius.  These  strains  can  than  be  used  to  measure  material  stretching,  petal 
displacement/curvature,  bending  work  for  the  final  deformed  geometry.  This  total  work  can 
than  be  loosely  related  to  the  energy  released  during  the  first  several  milliseconds  of  an 
explosion. 

Additionally,  a  series  of  successive  photographs  captures  the  formation  and  propagation 
of  radial  cracks  through  the  mesh.  In  doing  so,  one  can  estimate  the  incremental  strain  on  a 
given  mesh  element  and  estimate  Crack  Tip  Opening  Angle  (CTOA)  throughout  the 
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propagation  phase.  When  joined  with  a  known  force-displacement  history,  one  can  begin  to 
draw  conclusions  about  the  material’s  ability  to  resist  blast  damage. 

Numerical  Approach 

In  order  to  gain  insight  into  the  fracture  and  displacements  produced  expected  the  quasi¬ 
static  dishing  experiment,  a  numerical  model  was  created  in  ABAQUS.  An  axi-symmetric 
plate  model  was  statically  loaded  with  full-plate  and  half-radius  pressures  of  varying 
magnitude.  This  analysis  showed  the  (i)  mode  and  location  of  fracture,  (ii)  whether  fracture 
occurs  buy  thinning  or  necking,  and  (iii)  an  approximate  deformed  plate  shape  at  fracture. 


q  * 

$  T.~- 

OH]  IT  1  ITTT1 _ 5 

3 . t 

a  t  b 

R 

Figure  3:  ABAQUS  Numerical  Model  Representation 

Additionally,  a  moire  comprehensive  (and  considerably  more  accurate)  model  was 
developed  in  LS-DYNA  using  displacement  boundary  conditions  for  a  spherical  punch. 

Approximate  Analytical  Solutions 

Approximate  analytical  methods  have  been  previously  developed  for  each  of  the  quasi¬ 
static  phases.  Simonsen  (2000)  offers  an  approximate  closed-form  solution  (under  several 
assumptions)  for  plate  displacement  under  hemispherical  punch  loading.  With  that  solution, 
strains  and  dishing  work  can  be  computed  with  reasonable  accuracy.  These  results  are  then 
compared  with  those  obtain  experimentally  and  numerically. 
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Figure  4:  Generalized  Force  vs.  Displacement  History 

With  a  closed-form  solution  derived  originally  for  plate  cutting,  Wierzbicki  (1993)/(1999) 
developed  expressions  for  bending  and  membrane  work  done  during  the  circumferential 
cracking/ petalling  phase.  Combined,  the  various  methods  provide  a  reasonable  estimate  of 
the  total  work  done  in  deforming  the  plate  during  a  close-proximity  explosion. 

Still,  knowing  the  work  of  macro-scale  deformation  is  less  important  than  identifying  the 
overall  force-displacement  history.  The  deformed  plate  is  the  most  obvious  physical  clue  in 

i 

surmising  the  wave  front-plate  interaction  from  the  time  of  impart  through  to  the  petalling 
phase.  It  is  this  load-deformation  relationship,  combined  with  a  suitable  crack  propagation 
criterion,  which  can  most  closely  link  the  quasi-static  conditions  of  a  punch  test  to  the 
dynamic  conditions  present  in  an  explosive  blast. 


Approximate  Analysis 

Hemispherical  Dishing  and  Circumferential  Cracking  Phase 

As  mentioned  earlier,  the  dishing  analysis  is  based  on  work  done  by  Simonsen.  The 
theory  assumes  an  axi-symmetric  plate  is  dished  by  a  hemispherical  punch  of  radius  Rb.  A 
Cartesian  (w,  r)  coordinate  system  is  used  with  tpc  being  the  angle  from  the  center  of  the 
punch  to  the  outermost  contact  point  C.  The  total  punch  displacement  is  6  and  the  punch 
force  is  P.  Figure  5  illustrates  the  geometry  of  the  problem. 


Figure  5:  Spherical  Dishing  Geometry 

Several  assumptions  are  made  in  order  to  reach  a  closed-form  solution.  Plate  bending  is 
neglected,  plate  elements  are  considered  to  be  displaced  vertically,  and  the  material  is 
assumed  to  be  rigid-plastic.  Under  the  assumption  of  plane  strain,  the  generalized  radial 
membrane  force  is  given  by 


The  punch  force  (P)  at  point  C  can  be  expressed  as: 

P  =  2 7tN0Rb  sin2  y/(r) 

where  y/(r)  is  the  wrapping  angle  function 
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Given  the  relation  that: 
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The  following  differential  equation  can  be  surmised: 

dw  _  sin2^f 
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By  applying  the  boundary  condition  w{r)  =  0  at  r  =  R,  a  solution  for  w(r)  can  be  reached 
by  separation  of  variables: 
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for  rc  <  r  <  R 


(5) 


Consequently,  using  the  previously  stated  assumptions,  one  can  determine  the  vertical 
displacement  for  all  points  on  the  plate  and  the  force  required  to  deform  the  plate.  Of 
course,  at  some  point  the  plate  the  plate  will  reach  its  material  limits.  Strain  hardening 
characteristics  must  be  introduced  and  the  true  stress  strain  curve  is  assumed  to  obey  a 
power  law: 


To  predict  necking  localization  and  fracture,  the  maximum  on  the  load-displacement 


dp 

curve  is  identified  by  setting - 

dy/c 


=  0  .  One  can  then  solve  for  an  approximate  solution  for 


<pc  at  necking  failure.  The  final  expression  given  by  Simonsen  is: 


¥c_jm  *-957  +  .399h  (7) 

Using  this  relation,  additional  expressions  were  derived  to  determine  total  displacement  (Sfai]) 
and  total  work  (E)  up  to  the  point  of  necking  failure: 

6M,=\Aln*RMR*  (8) 
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Conical  Dishing 

In  some  cases,  a  cone  was  used  to  dish  a  flat  plate  with  a  pre-cut  starter  hole  in  order  to 
gain  insight  as  when  and  to  how  the  radial  cracks  form  and  propagate. 


Figure  6:  Conical  Dishing  Experiment 
The  problem’s  geometry  is  defined  in  Figure  7. 


Figure  7:  Conical  Dishing  Geometry 
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Again,  the  material  is  assumed  to  obey  a  true  stress  power  law. 
a  =  CJ" 


(10) 


Using  the  theory  of  moderately  large  deflections  in  thin  plates  and  following  Simonsen’s 
assumption: 


1  (  dwY 
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Consequently, 


The  generalized  membrane  stress  then  becomes: 


(12) 
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The  governing  equation  is  then: 
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Subject  to  the  boundary  conditions,  this  equation  has  the  general  solution: 

2» 

w  =  c,r2"+1  +  c2 


(15) 


The  constants  c,  and  c2  can  be  found  be  applying  the  boundary  condition 
r=rj.  The  expression  for  w  then  becomes: 
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where  £  is  the  radius  of  the  expanding  plastically  deformed  2one. 
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The  kinematic  boundary  condition  (K.B.C.)  is  then  applied  at  both  r  =  £  and  r  =  r,. 

v  r,  dw  d£  ,  n  , 

K  :  ^jT  +  ~dTW  =0  at  r=r>  (note  tv  =  -tan 0=  constant)  (17) 

where  T  is  time  and  r,  can  be  considered  a  time-like  parameter. 

Using  this  kinematic  condition,  £  can  now  be  expressed  in  terms  of  r,: 
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As  a  result,  (16)  becomes 
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The  vertical  punch  force,  P,  is  given  as  the  following: 

P  =  2?r,  JVtan  9  =  -?'|C‘>^tan  *£ — 

2" 


The  punch  displacement  is: 
£  =  +  r,  tan  #  =  2r,  tan  # 


(20) 


(21) 


Substituting  (20)  into  (21),  the  force-displacement  relation  is  developed. 


p_^C0/(tang)2wJ 

2" 

As  an  example,  consider  the  following  values: 

C0  =460  MPa  0  =30°  n  =  .22  /  =.045” 


P=48.5kN 


(22) 
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Force  (kN) 


The  analysis  provides  the  force-displacement  approximation  show  below  in  Figure  8. 


Approximate  Load  vs.  Displacement  Curve 

for  Varying  Thickness  Flat  Plates  using  120  °Cone 
80  - - 


0  10  20  30  1  40  50  61 


Displacement  (mm) 

Figure  8:  Analytical  Approximation  for  Conical  Dishing  with  Starter  Hole 


Radial  Cracking  and  Petalling 

The  theory  used  in  this  section  is  derived  in  detail  by  Wierzbicki  [1999],  He  proposes  that 
the  total  petalling  work  is  due  to  crack  propagation,  petal  bending,  and  membrane 
deformation.  Further,  these  quantities  are  shown  to  be  interdependent. 


Figure  9:  Theoretical  Petalling  Geometry 
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Due  to  physical  limitations  of  the  testing  apparatus,  fully  dished  plates  could  not  be 
pet  ailed  during  the  senes  of  experiments.  Consequently,  a  flat  plate  was  used  to  approximate 
the  deformation  and  fracture  conditions  seen  during  petalling.  Figure  6  shows  an  example 
of  a  six-petal  geometry  (npelal=6).  From  the  figure,  the  central  angle,  0,  is  given  as: 


^ petals 


(23) 


and  the  instantaneous  crack  length,  a,  is  related  to  the  petal  length,  1,  by 

/ 

a  = - 

cos# 


(24) 
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The  material  is  assumed  to  be  rigid-plastic  with  an  average  flow  stress,  oQ.  The  fully 


plastic  bending  moment  per  unit  length,  M0,  is  calculated  as: 


Ultimately,  the  rate  of  bending  work  in  a  single  petal  can  be  computed  as 

Wb  =  4M— /tan# 

P 

where, 

M  =  rjM0 


(25) 


(26) 


(27) 


M  is  the  amplified  bending  Moment,  which  accounts  for  the  increased  bending  resistance  of 
a  curved  plate.  In  the  flat  plate  approximation,  rf=  1.  The  per  petal  membrane  work  is  then 
shown  to  be  a  function  of  crack  tip  opening  displacement,  8t: 


^  _  3.84^/ V/3 /(sing)-473 
tcosQ 

where  p  is  the  instantaneous  bending  radius  of  the  flaps. 
Adding  (26)  and  (28)  to  get  total  work: 

w,=wb  +  wm 

then, 
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The  total  work  is  then  obtained  by  integrating  the  above  expression  with  respect  to  1: 
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Load  vs.  Displacement  for  .035"  Flat  Plate 


Figure  10:  Approximate  Theoretical  Load  Displacement  Curve  for  Petalled  Plate 

Punch  displacements,  the  measured  radius  of  curvature,  the  experimental  crack  tip 
opening  distance  were  all  used  in  estimating  punch  force.  Although  a  closed  form  expression 
for  P  vs.  1  was  not  developed  for  the  load  condition,  forces  were  estimated  using  changes  in 
total  work  over  small  increments  of  displacement. 
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(Detailed  computations  can  be  reviewed  in  Appendix  G.)  In  this  case  a  piece-wise  force- 
displacement  curve  for  each  petalled  plate  was  generated  using  Wierzbicki’s  expressions  and 
experimentally  measured  inputs.  An  example  of  a  generated  curve  is  shown  in  Figure  10. 


Strain  Field  Development 


!• 


Figure  12:  Incremental  Strain  Field  Rotation 

As  shown  in  Figure  12  above,  the  strain  field  can  then  be  rotated  into  alignment  with  the 
direction  of  crack  growth  by  multiplying  by  the  appropriate  direction  cosine.  Although 


principal  stresses  are  not  know,  the  stress  ratios  (specifically,  which  governs  necking) 

a22 

can  be  determined  from  en  and  e22  using  the  Von  Mises  yield  condition  and  the  associated 
flow  rule. 
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Divide  (36)  by  (37)  and  solve  for  —  in  terms  of  to  obtain: 
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Numerical  Modeling 

Preliminary  Dishing  Model  in  ABAQUS 

ABAQUS  was  used  to  model  the  preliminary  plate’s  loading  condition.  For  convenience, 
a  pressure  loading  condition  (vice  a  displacement  boundary  condition)  was  selected  as  a  first 
approximation.  In  this  model,  a  uniform  pressure  was  first  distributed  over  the  entire  plate. 
This  loading  geometry  resulted  in  dishing  without  reverse  curvature  and  a  necking  fracture  at 
the  plate  s  clamped  boundary.  The  uniform  pressure,  full  radius  loading  analysis  was  not 
pursued  further,  since  the  load  condition  result  in  a  realistic  deformation  mode. 


Radius  (mm) 


Figure  13:  Displacement  vs.  Radius  for  Half-Radius  Pressure  Loading 


Figures  14  and  15  show  the  displacements  for  the  half  radius  pressure  distribution.  It 
was  found  that  using  a  pressure  distribution  over  half  the  plate’s  radius  more  closely 
modeled  the  conditions  of  hemispherical  punch  dishing. 


Figure  14:  Sectional  View  of  Deformed  Plate  at  240  MPa  Pressure  Amplitude 


Figure  15:  View  of  Deformed  Plate  at  240  MPa  Pressure  Amplitude 


LS-DYNA  Model 


Towards  the  end  of  testing,  a  Finite  Element  model  was  created  in  LS-DYNA  using  the 
displacement  boundary  conditions  of  a  rigid  sphere,  shown  in  Figure  16.  Though  the 
revised  model  moved  further  from  the  conditions  of  gas  pressure  loading  (as  in  the 
preliminary  simulation),  it  more  accurately  represented  the  quasi-static  experiment.  Three 
conditions  were  modeled:  (1)  the  75  mm  radius  sphere/  .055”  plate,  (2)  50  mm  radius 
sphere/,045”  plate,  and  (3)  the  20  mm  radius  sphere/  .035”  plate.  Computational  results 
were  then  compared  with  both  the  approximate  analytical  predictions  (from  Simonsen)  and 
the  experimental  values.  Data  from  this  analysis  can  be  found  in  the  “Results  and 
Discussion”  section  of  this  report.  A  sample  LS-DYNA  input  file  for  theR=75mm/.055” 
plate  test  case  can  be  found  in  Appendix  H. 
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Figure  16:  LS-DYNA  Simulation  with  Rigid  Body  Sphere  Boundary  Condition 


Design  and  Testing 


Apparatus  Design 

A  specialized  fixture  was  fabricated  for  indentation  testing  of  thin  square  sheets. 
Designed  to  work  using  a  200  KN  load  cell  in  a  Universal  Testing  Machine,  the  fixture 
(shown  below  in  Figure  17)  can  safely  accept  a  central  point  load  (P)  up  to  the  200  KN 
machine  limit.  Edge  fixity  is  achieved  through  the  use  of  serrated  bolting  ring  mating 
surfaces.  Detailed  size  and  material  specifications  for  the  fixture  are  shown  in  Figure  18. 
Further  design  explanations  can  be  found  in  Appendix  D. 


Figure  17:  Notional  Test  Fixture  and  Completed  Design 

To  create  the  dishing  effect,  three  hemispherical  indenters  were  used  with  radii  of  20  mm, 
50  mm,  and  75  mm.  Indenter  speed  (V)  was  set  at  a  baseline  of  10  mm/min.  Indentation 
was  continued  through  circumferential  cracking  until  the  instance  at  which  radial  fracture 
began.  At  that  point,  the  dishing  portion  of  the  test  was  stopped.  A  120°  conical  indenter 
was  then  used  to  simulate  petalling  from  the  flat  plate  condition  for  the  same  thickness. 
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Specimen  Design 


Figure  19:  Test  Plate 

Figure  19  shows  a  representative  test  plate.  Three  plate  thickness  were  tested  (.035”, 
.045”,  and  .055”).  The  plate  material  has  an  average  uniaxial  tensile  yield  strength  (ay)  of 
29,000  psi  (200MPa),  an  average  ultimate  tensile  strength  (aj  of  46,000  psi  (317  MPa),  and  a 
Rockwell  Superficial  Hardness  of  72.  Detailed  matenal  characteristics  can  be  reviewed  in 
Appendix  E.  The  uncoated  side  of  the  plate  is  scribed  with  a  2.5  mm  strain  mesh  (slightly 
larger  than  2  plate  thicknesses).  During  deformation,  the  horizontal  and  vertical  local 
incremental  strains  are  directly  recorded  using  digital  photography. 
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Testing 


Table  1  shows  the  test  matrix  of  the  series  of  experiments  performed  with  the  Universal 
Testing  Machine.  Clearly,  there  were  many  testing  variables  to  consider.  With  a  finite 
number  of  test  specimens,  however,  it  was  a  goal  of  the  study  to  develop  a  testing 
methodology  and  identify  the  most  significant  factors. 

Tests  1  through  3  were  spherical  indentation  tests,  which  were  used  to  verify  Simonsen’s 
predictions.  Tests  4  through  7  were  used  to  determine  the  best  approach  given  the 
limitations  on  time  and  equipment.  In  Test  4,  a  conical  punch  was  used  to  indent  a  plate 
with  a  preformed  hole.  Later,  Tests  14  andl5  were  used  to  observe  to  effect  of  hole  size  in 
this  same  loading  condition.  Test  5  was  the  first  experiment  in  which  the  starter  hole  was 
pre-cracked  in  order  to  promote  the  petalling  behavior.  Tests  8,  9,  and  10  also  used  pre¬ 
cracked  holes  to  examine  the  effect  of  plate  thickness.  Tests  11  and  12  investigated  the 
effect  of  increasing  the  number  of  pre-cracks,  while  Test  13  was  used  to  observe  the  effect 
of  increasing  pre-crack  length.  Tests  4, 14  and  15  confirmed  the  effect  of  starter  hole  size. 
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Results  and  Discussion 


Spherical  Indentation 

Figure  19  shows  the  experimental  and  approximate  analytical  results  for  spherical 
indentation  testing.  (Note,  the  analytical  approximation  is  valid  only  up  to  Simon  sen's  failure 
prediction.  Points  beyond  the  maximum  load  are  artifacts  of  the  computation.) 


Spherical  Indentation  Results 


Cross  Head  Displacement  (mm) 


Figure  20:  Approximate  and  Experimental  Results  for  Spherical  Indentation 

Figures  21  through  23  on  the  following  page,  show  results  of  the  LS-DYNA  numerical 
simulations  for  the  same  three  load  cases. 


Spherical  Indentation  Results  for  R=20mm  Sphere  /  .035"  Plate 
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Figure  21:  LS-DYNA  Results  for  R=20  mm  punch  /.035”  Plate 


Spherical  Indentation  Results  for  R=50  mm  Sphere  /  .045"  Plate 
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Figure  22:  LS-DYNA  Results  for  R=50  mm  punch  /.045”  Plate 


Spherical  (dentation  Results  for  R-75  mm  Sphere  /  .055"  Plate 
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Figure  23:  LS-DYNA  Results  for  R=75  mm  punch  /.055”  Plate 


m 


As  in  Simonsen’s  experiment,  the  experiments  correlated  to  the  analytical  approximation 
within  5%.  The  preliminary  LS-DYNA  model  (pictured  in  Figure  24)  also  resulted  in  good 
agreement,  although  the  model  requires  further  refinement  for  follow-on  tests.  This 
portion  of  testing  validated  the  testing  method  and  clearly  proved  that  the  new  specimen  test 
fixture  functioned  as  designed. 

During  spherical  indentation  testing,  it  became  clear  that  the  specimen  material  was  too 
ductile  to  carry  out  all  tests  as  planned.  Plastic  deformations  were  quite  large  during 
dishing,  and  the  shallow  conical  punches  did  not  provide  sufficient  depth  of  stroke  to 
continue  petalling.  Furthermore,  the  dishing  failure  always  resulted  in  incomplete 
circumferential  cracking  (as  in  Figure  24  below).  As  a  result  of  these  two  obstacles,  the 
dishing  and  petalling  phases  were  modeled  in  two  entirely  separate  experiments. 
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150  mm  Diameter/.055”  Plate 


40  mm  Diameter/.035”  Plate 


Figure  24:  Experimental  and  LS-DYNA  Spherical  Dishing  Deformations 
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Conical  Indentation  without  Pre-cracks 


Figure  25  shows  the  force-displacement  relation  observed  during  conical  dishing.  During 
this  experiment,  a  constant  plate  thickness  was  tested  while  the  starter  hole  was  varied  in 
size.  The  dashed  line  represents  the  analytical  approximation  derived  earlier.  The  analytical 
approximation  suggests  that  the  force-displacement  curve  should  be  independent  of  hole 
size.  This  behavior  was  observed  up  to  a  displacement  of  about  20  mm.  After  that  point, 
the  3  experiments  diverge.  Certainly  the  smaller  hole  size  induced  much  earlier  radial 
cracking  due  to  much  higher  hoop  stresses. 

It  is  difficult  to  draw  immediate  conclusions,  since  each  case’s  geometry,  and 
consequently  the  material’s  load  history,  cannot  be  directly  correlated.  The  analytical 
approximation,  however,  serves  as  a  broad  estimate  of  the  material’s  general  force- 
displacement  behavior,  but  falls  short  in  accounting  for  the  problem’s  changing  geometry. 


.035"  Plate  Dished  with  120°  Cone 

(10,  22,  and  45  mm  Hole,  No  Precracks) 


Figure  25:  Conical  Indentation  Force-Displacement  Curves  (with  Starter  Hole) 


Conical  Indentation  With  Pre-cracks  (Petalling) 


Load  vs.  Displacement  for  Varying  Thickness  Flat  Plate 
120°  Cone,  Four  3-mm  Pre-cracks,  9.5  mm  Hole 


Initial  Condition:  9.5  mm  Hole  with  4  Pre-cracks 
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Figure  26:  Conical  Indentation  Force-Displacement  Curves  (with  Pre-cracks) 


Figure  26  shows  the  results  of  the  plate  petalling  experiment,  which  was  conducted  with  a 
9.5  mm  starter  hole  and  four,  3  mm  pre-cracks.  Figure  27  shows  how  the  petalling  test 
progressed. 


10  mm 


25  mm 


55  mm  (completed  test) 


Figure  27:  Petalling  Test  Punch  Displacement  Progression 


The  petalling  test  showed  that  the  effect  of  plate  thickness  was  approximately  linear.  In 
each  experiment  during  the  first  10  mm  of  displacement,  the  cone  settled  into  the  starter 
hole  and  propagated  the  pre-cracks.  After  that  point,  the  petals  began  to  form.  From  the 
pictures,  one  can  see  that  very  little  curvature  developed.  Again,  this  was  due  to  the  shallow 
angle  of  the  conical  punches.  Only  55  mm  of  punch  displacement  could  be  observed  before 
the  punch  began  to  interfere  with  the  fixture’s  boundary  conditions,  at  which  point  the  test 
had  to  be  stopped. 

Although,  there  appears  to  be  disagreement  between  the  analytical  approximation  and  the 
observed  force-displacement  curves,  the  general  trend  and  level  of  force  correlate  reasonably 
well.  It  is  apparent  that  the  artificiality  of  pre-cracking  has  some  impart  on  the  experiment. 
In  the  two  experiments  where  natural  radial  cracks  formed  and  petalling  could  be  observed 
(Figure  28),  the  force-displacement  curve  closely  resembled  the  behavior  derived  out  of 
Wierzbicki’s  expressions  Eq.(32). 


Displacement  (mm) 


Pre-Dished  and  Cut  Specimen  Conical  Dishing  with  Starter  Hole 

Figure  28:  Force-Displacement  Observations  During  Petalling 


It  is  important  to  keep  in  mind  that  the  dishing  force  is  five  to  six  times  higher  than  the 
less  critical  petalling  force.  With  the  perspective  that  this  approach  is  simply  another 
available  tool,  this  petalling  approximation  could  certainly  provide  useful  first-order 
estimates  to  a  designer  or  post-explosion  investigator. 

Effect  of  Pre-crack  Length 

Figure  29  shows  the  effect  of  pre-crack  length.  Variation  in  the  length  of  the  starter 
crack  produced  little  effect,  other  than  lengthening  the  amount  of  time  required  for  the  cone 
to  settle  into  full  contact  with  the  plate.  The  result  was  a  small  shift  in  the  force  - 
displacement  curve.  As  displacement  increased  the  effect  was  diminished.  This  behavior 
was  expected,  but  required  verification.  The  length  of  pre-crack,  however,  did  affect  the 
number  of  petals  formed.  The  3-mm  pre-crack  yielded  four  petals,  while  the  10-mm  pre¬ 
crack  yielded  only  four.  This  behavior  agrees  with  Wierzbicki’s  prediction  than  some  energy 
absorption  minimum  exists,  tending  to  promote  the  formation  of  three  to  five  petals. 
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Figure  29:  Effect  of  Pre-Crack  Length  on  Force-Displacement  Curve 


Effect  of  Number  of  Pre-cracks 

As  shown  in  Figure  30,  the  number  of  pre-cracks  had  no  significant  effect  of  the  overall 
force-displacement  curve.  Only  four  petals  developed  in  each  case,  however.  In  both  tests, 
all  pre-cracks  tended  to  start,  but  only  four  continued  to  grow.  As  in  the  previous 
examination  of  pre-crack  length,  this  result  also  suggests  that  a  petalling  energy  minimum 
exists,  which  is  independent  of  the  number  of  initial  radial  necks  or  cracks. 


Effect  of  Number  of  Precracks 
3  mm  Pre-cracks,  .035"  Plate,  9.5  mm  Starter  Hole 


Cross  Head  Displacement  (mm) 


Figure  30:  Effect  of  Number  of  Pre-Cracks  on  Force-Displacement  Curve 
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Crack  Displacement  and  Strain  Field  Examination 


Figure  31  shows  the  increment  of  crack  growth  used  while  mapping  the  strain  field 
around  the  crack  tip.  The  crack  was  allowed  to  grow  approximately  2.5  mm  between  framed. 
The  white  arrow  gives  a  fixed  point  from  which  to  visually  reference  the  crack’s  advance. 

The  mapping  grid  was  laid  out  according  to  the  method  discussed  in  an  earlier  section. 


Frame  1  Frame  2 

Figure  31:  Crack  Advance  Increment 


The  displacement  filed  vectors  are  plotted  relative  to  the  crack  tip.  Examination  of  the 
displacement  field  shows  that,  at  this  instant  in  time,  the  membrane  is  rapidly  stretching  and 
thinning  in  the  region  in  front  of  the  crack. 

Biaxial  incremental  strain  values  were  computed  for  a  small  rectangular  patch  around  the 

crack  tip.  Similarly,  the  principal  stress  ratio,  — ,  was  computed  from  Eq.(38)  for  each  of. 

<J22 

the  grid  points,  helping  to  identify  regions  where  necking  becomes  critical.  Figure  32  shows 
the  measured  displacement  field  in  the  postage  stamp  size  area  surrounding  the  crack  tip. 
Figure  33  shows  the  measured  incremental  strain  fields  and  Figure  34  shows  the  computed 


field.  A  visual  inspection  of  the  displacement  and  strain  fields  suggests  that  necking 

fj22 

occurs  primarily  directly  ahead  of  the  crack.  This  behavior  was  both  expected  and  observed 
during  the  experiment. 
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Figure  32:  Crack  Advance  Incremental  Displacement  Field 
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Figure  33:  Resultant  Strain  Increment  Fields 
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Conclusions  and  Recommendations 


Conclusions 

This  study  has  investigated  using  simple,  inexpensive  quasi-static  tearing  tests  to  gain 
insight  into  the  complex  dynamic  problem  of  explosions.  It  has  been  shown  that  these  tests 
can  successfully  reproduce  plate  deformations  similar  in  nature  to  those  found  after  large- 
scale  explosions.  It  was  also  shown  that  the  forces,  displacements,  and  energy  dissipation  in 
the  plate-dishing  phase  can  be  accurately  approximated  by  Simonsen’s  approach  (within  5%). 
Although  this  series  of  experiments  was  not  able  to  generate  spiraled  petals,  Wierzbicki’s 
calculation  of  petalling  work  appeared  to  provide  reasonably  good  (within  20%)  estimates  of 
energy  dissipation  during  the  second  phase  of  deformation.  The  conical  punch  test, 
however,  requires  a  revised  mechanical  design  in  order  to  better  model  the  flat  plate  petalling 
and  post-dish  petalling  cases. 

An  approach  to  developing  fracture  criterion  was  also  presented  by  measuring 
incremental  strains  around  the  crack  tip.  Subsequently,  the  local  stress  ratios  governing 
material  necking  were  computed.  With  that  information  in  hand,  crack  growth  can  be 
predicted  in  the  next  increment.  This  approach  could  result  in  a  revised  crack  growth 
criterion,  which  might  then  be  used  to  improve  existing  hydrodynamically  coupled  finite 
element  simulations. 
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Recommendations 


•  The  punch  tests  should  explore  a  wider  range  of  penetrator  geometries.  Punch  variation 
may  introduce  more  curvature  into  the  resulting  petals,  helping  to  capture  more  of  the 
dynamic  effects.  Pyramid  or  volute  shapes  might  be  considered.  Certainly,  narrower 
cone  angles  should  be  tried  in  order  to  increase  the  punch  displacement  and  extend 
observation  time. 

•  Future  analytical  work  should  attempt  to  more  closely  link  the  quasi-static  test  results 
with  characteristics  of  the  dynamic  problem.  To  do  so  may  require  dynamic  testing, 
which  could  be  accomplished  at  government  blast  testing  facilities. 

I 

•  Future  numerical  simulations  should  apply  new  punch  geometries  and  include  new  plate 
structures  (sandwiches,  foam,  etc... )  which  might  potentially  be  considered  for  ship  hull 
reinforcement.  Future  numerical  simulations  should  also  attempt  to  more  completely 
capture  the  dynamic  aspects  of  an  explosion. 

•  The  potential  improvement  in  crack  growth  criterion  could  result  in  a  tremendous 
advance  m  the  understanding  of  crack  propagation  due  to  intense  blast  pressures.  This 
avenue  should  be  explored  more  fully  in  order  to  update  current  prediction  codes. 

•  Other  areas  of  application  for  this  research  should  be  investigated.  Certainty  pipeline 
explosions  (as  in  the  natural  gas  industry),  internal  aircraft  explosions,  and  pressure 
vessel  explosions  are  just  of  few  items  for  consideration 
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Appendices 


Appendix  A:  Introductory  Material 


Since  the  inception  of  the  warship,  designers  have  sought  means  to  make  combatant 
vessels  more  damage  resistant.  Only  since  World  War  II,  however,  has  there  been  significant 
effort  to  compile  and  analyze  detailed  reports  of  vessel  battle  damage.  In  the  subsequent 
post  war  years,  much  analysis  was  conducted  in  the  areas  of  hull  armor,  resistance  to  mine 
attack,  and  torpedo  side-protection.  As  might  be  expected,  the  systems  that  were  ultimately 
developed  to  enhance  survivability,  incurred  substantial  weight,  space,  and  performance 
penalties.  As  the  years  have  passed,  the  damage  database  has  deteriorated.  Since  the  late 
1940’s,  the  U.S  Navy  has  relied  largely  on  full-scale  tests  or  “SINKEX’s”  to  demonstrate 
battle  damage  effects.  There  have  been  only  a  handful  of  wartime  incidents  demonstrating 
the  realistic  damage  absorbing  and  recovery  capabilities  of  modem  warships.  In  that  time, 
threat  weapons  (mines,  missiles,  and  torpedoes)  have  been  improved  markedly.  In  contrast 
to  their  World  War  II  predecessors,  state  of  the  art  weapons  are  extremely  likely  carry  out  a 
successful  attack  and  inflict  extraordinary  damage  with  a  high  degree  of  accuracy. 

As  seen  more  recently  in  the  case  of  the  U.S.S.  Cole,  tremendous  hull  damage  can  be 
suffered  during  unconventional  close  aboard  explosions.  Figure  A1  shows  the  20-ft  by  40-ft 
amidships  hole  torn  in  the  ship’s  port  side  by  a  waterline  explosion. 

Inspection  of  the  picture  shows  that  the  explosion  created  a  spherical  bulge,  or  dishing,  in 
the  ship’  hull,  prior  to  tearing  into  numerous  petals.  Though  officially  unconfirmed,  it  has 
been  suggested  that  the  attack  was  carried  out  using  between  400  and  700  lbs.  of  C-4  military 
explosive.  Unofficial  standoff  distance  estimates  range  from  on-contact  to  10  feet. 


Figure  Al:  USS  Cole  Port  Side  Damage 


In  this  instance,  a  majority  of  the  damage  occurred  below  the  waterline.  Figure  2  clearly 
shows  the  above-waterline  tearing.  From  Figure  A2,  the  characteristic  petal  formations  car 
be  identified,  as  well  as  the  upper  deck  hard  point  at  which  crack  propagation  is  arrested. 
Hull  panels  at  this  location  on  a  typical  naval  combatant  ship  are  usually  constructed  of 
either  mild  steel  (MTS-45)  or  HY-80,  and  they  can  vary  between  .375”  and  .5”  in  thickness 


m 


Appendix  B:  Additional  Background 

Holing  prediction  for  metal  plate  has  been  a  topic  of  study  since  the  early  1900’s.  The 
research  is  rooted  in  a  1912  study  done  by  Bertram  Hopkinson,  which  examined  the 
resistance  of  naval  plating  against  artilleiy  shell  penetration.  In  that  case,  the  deformed 


Figure  Bl:  Hopkinson's  1912  Armor  Plate 

plating,  shown  in  Figure  Bl,  developed  a  coronet-shaped  pattern  with  radial  cracks  and 
necking  at  the  petal  edges.  Since  that  time  there  have  been  countless  experiments 
investigating  armor  effectiveness.  Until  the  past  two  decades,  the  vast  majority  of  plate 
holing  and  deformation  studies  have  focused  on  the  impact  and  associated  damage  of 
ballistic  particles.  The  most  notable  and  comprehensive  exception  is  the  1940’s  work  of  Sir 
G.I  Taylor  (1948)  who  examined  the  nature  of  submerged  blast  waves  and  their  effects  on 
thin  plates.  Taylor’s  work  in  the  area  of  submerged  explosions  was  continued  by  Robert 
Cole  (1948)  in  the  text,  Underwater  Explosions.  Together,  these  two  researchers  provide 
much  of  the  technical  underpinning  for  current  blast  damage  prediction  methods. 
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Explosive  deformation  and  holing  studies  in  naval  vessels,  largely  classified,  have  been 
empirical  derivations  based  on  years  of  accrued  data.  An  undisclosed  solution  developed  by 
ONR,  predicts  a  minimum  and  maximum  expected  hole  radius  in  naval  panels  for  a  given 
set  of  panel  parameters.  This  engineering  tool  was  the  result  of  hundreds  of  live  fire  tests, 
and  is  one  of  the  most  commonly  used  methods  used  by  U.S.  Navy  designers  in  estimating 
expected  battle  damage  due  to  close  proximity  explosions.  This  solution  suggests  the 
following  general  relationship: 


Stiffened  Panel 


Figure  B2:  Current  Panel  Damage  Formulation 


^Min  —  —  ^Max  =  /(G,  T,  M) 

R),  =  Resultant  Hole  Size 

R™,  =  Minimum  Predicted  Hole  Size 

R^  =  Maximum  Predicted  Hole  Size 

G  =  Panel  Geometry 

T  =  Plate  Thickness 

M  =  Material  Properties 
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Other  problems,  however,  lend  themselves  to  the  study  of  panels  under  explosive  loading. 
Wierzbicki  and  Thomas  (1993),  by  examining  plate  cutting  behavior,  were  able  to  related  the 
geometiy  of  the  deformed  plating  geometry  to  the  mechanical  and  frictional  work  done  in  a 
vessel  grounding,  thereby  predicting  damage.  Later,  Wierzbicki  (1996)  and  Nurick  (1996), 
investigated  the  response  of  clamped  thin  sheets  when  subjected  to  on-contact  charges. 
Wierzbicki  (1999)  subsequently  proposed  that  the  kinematics  of  the  cutting  process,  shown 
in  Figure  B3,  were  very  similar  to  both  those  of  the  explosive  petalling  and  ballistic 
penetrator  problem. 


Figure  B3:  Cutting  and  Petalling  Geometric  Similarities 

Along  similar  lines,  Atkins  (1998)  observed  the  necking  and  radial  crack  formation  in 
ductile  materials  when  perforated  by  both  spherical  and  conical  penetrators.  Using  a 
hydraulic  bulger,  he  later  analyzed  thin  sheet  necking  around  biaxially  loaded  holes  (1999). 
Simonsen  (2000),  built  upon  this  work,  using  spherical  indenters  to  analyze  force- 
displacement  relationship  in  clamped  mild  steel  plates.  Again,  these  studies  were  targeted  at 
making  damage  prediction  for  vessel  collisions  and  ballistic  penetrators.  Nevertheless,  these 


experiments  closely  model  thin  sheet  behavior  during  the  dishing,  disking,  and  pre-petalling 
phases  of  an  explosion. 

To  date,  there  has  been  little  effort  to  correlate  the  results  of  quasi-static  indenter 
experiments  to  the  damage  observed  both  in  live  fire  explosive  tests  and  in  combat.  With 
that  goal  in  mind,  this  research  is  deemed  both  relevant  and  unique. 


Appendix  C:  Problem  Statement  Details 


Currently,  U.S.  Navy  surface  ships  are  built  using  a  damaged  length  design  standard 
equivalent  to  15%  Length  Between  Perpendiculars  (LBP).  This  rule  for  damaged  length  is 
based  on  both  past  operating  experience  and  on  the  damage  expected  from  a  given  threat 
weapon.  As  illustrated  in  Figure  Cl,  accurately  determining  the  opening  size  is  crucial  in 
determining  compartment  layout,  equipment  placement,  and  bulkhead  spacing.  In  a  poorly 
designed  500-ft  ship,  for  instance,  a  15%  damage  length  (75ft)  could  result  in  total  loss  of 
propulsion.  By  the  same  token,  over  designing  a  ship  using  excessive  damage  length  results 
in  increased  structural  weight,  higher  cost,  and  decreased  performance. 


Figure  Cl:  Surface  Ship  15%  LBP  Damage  Length 

Damage  prediction  methods  can  be  extraordinarily  complex.  The  most  modem 
techniques  use  machine  codes,  which  couple  explosive  mechanics  with  the  concurrent 
hydrodynamic  effects.  These  codes,  however,  are  both  expensive  and  time  consuming. 
Furthermore,  the  accuracy  of  their  output  is  not  yet  fully  validated.  Other  than  scaled  live 


fire  tests,  no  reliable  prediction  tool  exists  which  can  quickly  and  accurately  estimate  damage. 
An  analytical  method  which  uses  inputs  such  as  charge  size,  standoff  distance,  material  could 
prove  useful  both  to  ship  and  weapon  designers. 


Appendix  D:  Fixture  Design  Supplemental 


A  significant  amount  of  design  work  was  required  to  ensure  that  the  mechanism  (1)  did 
not  fail  under  the  expected  test  loads  and  (2)  was  adaptable  for  use  in  the  Universal  Testing 
Machine  (UTM)  belonging  to  the  MIT  Impact  and  Crashworthiness  Laboratory. 

As  shown  in  the  drawings,  the  fixture  horizontally  holds  a  300  mm  x  300  mm  square  thin 
sheet  while  an  indenter  or  punch  is  lowered.  The  clamped  dimensions,  shown  in  Figure  Dl, 
are  220  mm  x  220  mm.  The  specimens  maximum  dimensions  were  selected  based  on 
several  criteria.  First,  the  specimen  needed  to  be  large  enough  to  easily  observe  the  resultant 
deformations,  yet  small  enough  to  fit  within  the  UTM.  Further,  it  was  desirable  to  have  a 
portable,  light-weight  test  fixture  (in  this  case  75  lbs.  was  designated  as  an  upper  weight 
limit).  Lastly,  steel  sheet  stock  is  normally  manufactured  in  4  ft  x  8  ft  sheets.  A  300  mm 
square  sample  allows  32  specimens  to  be  taken  from  a  single  sheet  with  minimal  waste. 


220  mm  Clamped  Width 


I 


1.5  mm  Fillet  Radius 


300  mm  Specimen  Width 
(1.15  mm  Specimen  Thickness) 


10  mm 


25  mm 


Figure  Dl:  Clamped  Specimen  Cross  Section 


Since  the  test  specimen  represents  a  single  clamped  panel,  the  specimen  needed  to  be 
geometrically  similar  to  existing  naval  panels.  Frame  spacing  within  a  mid-sized  combatant 
is  generally  near  8  ft,  with  the  panels  being  nearly  square  (or  slightly  rectangular).  At 
thicknesses  of  1/2”  to  3/8”,  a  panel  width-to-thickness  ratio  (w/t)  of  200  can  be  reasonably 
assumed.  In  order  to  maintain  this  ratio,  three  specimen  thicknesses  (.89mm,  1.14  mm,  and 
1.40  mm  )  were  selected  from  available  mild  steel  stock. 

The  body  of  the  fixture  is  machined  from  mild  steel  (45ksi  Yield  Strength).  Member 
thicknesses  were  selected  to  prevent  yield  during  normal  testing  and  to  minimize  weight. 
Note  that  the  fixture  narrows  to  6”  at  its  base  to  transfer  force  directly  to  the  load  bearing 
beam  in  the  UIM.  Thin  lateral  bracing  panels  were  used  to  prevent  fixture  racking  during 
testing.  The  fixture  mounts  to  the  UTM  through  a  series  of  threaded  holes  in  its  baseplate. 

The  upper  and  lower  bolting  rings  are  used  to  hold  the  specimen  in  place.  These  rings  are 
machined  from  A2  Air  Hardened  Steel  (108  ksi  Yield  Strength).  This  grade  of  steel  was 
selected  to  satisfy  fixity  requirements.  Firm  edge  fixity  is  required  in  order  to  make 
boundaiy  condition  assumptions.  Past  experiment  have  shown  that  simple  bolting  flanges 
do  not  provide  sufficient  clamping  restraint.  Typically  a  drawbead  can  be  used  to  overcome 
this  problem.  When  testing  thick  or  high  strength  materials,  however,  it  becomes  difficult  to 

j 

form  the  specimen  around  a  drawbead.  For  this  reason  (and  since  the  fixture  is  intended  to 
test  a  wide  range  of  materials  and  thickness),  a  serrated  mating  surface  was  used  to  provide 
fixity. 

In  order  to  ensure  sufficient  tooth  hardness  and  clamping  force,  an  estimate  of  material 
strength  was  made  using  Slip  Line  Fracture  Mechanics.  The  analysis  showed  that  the  mating 
surfaces  must  have  a  yield  strength  at  least  50%  higher  than  that  of  the  test  specimen.  The 
estimates  are  summarized  in  Figures  D2  and  D3. 


Test  Material 

.  X 


,45  Degrees 


mm 


Figure  D2:  Mohr  Circle  For  Tooth  Element  Tooth 


■crrf=-tf-k=-2k((l+^) 


tfrW  =  -2kw(l+f) 


d  =  -k-2k 


Figure  D3:  Mohr  Circle  For  45-Degree  Rotation  Tooth  Element 


=  ^n, 


M,+f>  =  Mi+f) 

k  ]  +  ~ 

T~  =  “■  =  1  -44  (Ratio  of  Base  Plate  Yield  to  Test  Material  Yield) 


Estimate  of  Frame  Sizes  and  Stresses  Induced  During 
Indentation  Testing 


•  Cross  Beam  Dimensions: 

h  =  25.4mm  b  :=  50- mm  L:=  300mm 

Expected  Load  and  Material  Characteristics: 

P  :=  10500CN  P  =  2.36  x  104  lbf 

Y  :=  1 08000 psi  Y  =  7.446x  108Pa 


Assuming  a  Rigid  Test  Specimen  and  a 

Simply  Supported  Beam  with  a  Concentated  Central  Load: 


M  := 


L  P 
2  4 


abeam  ‘  j  abeam  _  ^.324x  10  Pa 

Safety  Factor: 

Y 

- =  1.017 

abeam 

...or  distibuting  load  along  ther  beam,  gives  the  same  result  for  the  maximum 
developed  stress: 

i 

P  ^  1  T2 

w:= —  M  :=  —  •  vv-L 

4L  2 

°beam‘=^’Y  cbeam~^'^^x  ^  Pa 

Safety  Factor: 

Y 

- =  1.017 

abeam 
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Using  the  exact  solution  for  sinusoidally  loaded,  simply  supported  plate: 

FD;=-7  M:=FDi 

71 

CTbeam  :=  M  •  — 

Safety  Factor: 

Y 

- =  2.509 

abeam 

Compressive  Load  on  Legs: 

Yjeg  :=  45000  psi 

w,eg  :=  49.3  mm  Aleg  :=  2(2t|eg- wIeg) 

<T|eg  =  2.662  x  107  Pa 


CTbeam  —  2-968  x  10  Pa 
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Appendix  E:  Uniaxial  Testing 


The  mild  steel  plate  material  was  tested  Uni. -axially  in  tension  according  to  ASTM 
standards.  The  .035”  and  .045”  were  tested  in  1999  in  concert  with  a  separate  project.  The 
.055”  thickness  was  tested  specifically  for  this  research.  All  specimens  were  tested  in  the  X 
and  Y  direction  to  verify  that  the  material  behaved  isotropically. 


Load  vs.  Time  for  .035"  X  Sample  2 


125  175  225  275  325  375  425 


Time  (Sec) 

Figure  El:  Representative  Load  vs.  Time  Curve 


Load  Change  vs  Time  .035  X  Sample  2 


Figure  E2:  Representative  Load  Change  vs.  Time  Curve 


During  the  uni-axial  testing,  all  of  the  thin  strip  specimens  displayed  stable  fracture, 
consequently  the  precise  load  at  fracture  was  not  readily  apparent.  By  plotting  the  change  in 
load  vs.  time  for  each  specimen,  a  “knuckle”  in  the  load  vs.  Time  plot  was  clearly  identified, 
which  indicated  the  onset  of  fracture.  Figures  El  and  E2  illustrate  this  approach.  Figures 
E3  and  E4,  respectively,  show  a  representative  uni-axial  specimen  and  a  close-up  view  of  the 
fracture. 


The  final  dimensions  of  each  specimen’s  cross  section  were  measured  in  order  to 
determine  the  true  stress  and  strain  of  the  specimen  at  fracture.  This  data  point  was  vital  in 
constructing  the  true  stress-strain  curve  for  the  finite  element  simulation.  Due  to  necking, . 
the  specimen’s  highest  true  fracture  stress  and  strain  were  achieved  in  the  middle,  where  the 
fracture  initiated.  Figure  E5  shows  the  approximate  final  geometry  for  the  necked 
specimens. 


Figure  E5:  Approximate  Necked  Specimen  Geometry 


C* 

Reduction  in  Squared  Thickness  (RST)  and  true  fracture  strain  for  the  thin  strips  were 
computed  in  accordance  with  McClintock  and  Zheng  (1993), 


RST  =  \-{tflt0f 
ss  =  -ln(\- RST) 


(El) 

(E2) 
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m 


dtf 


Figure  E6:  Representative  True  Stress  vs.  Strain  Curve  for  Uniaxial  Test 


Some  formulation  was  required  to  compute  onilM[  (o^J  : 


£e)~r 

as 


(E3) 


P fracture  (E4) 

Where  substituting  (E3)  into  (E4)  yields, 

^ 'fracture  ~  ^m^ni  4"  (®"m  —  (^m  ~  ^e)~l  ) 


m 

(E5) 


and  solving  (E5)  for  om 

rr  _ _ _ as 


A,  +  A 


(E6) 


Table  El:  Edge  and  Middle  Thickness  at  Fracture  for  Uniaxial  Tests 
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m 


Table  E2:  Right  and  Left  Edge  Fracture  Strain  for  Uniaxial  Tests 


Sample 


.035  X#1 


.035  X  #2 


Right  Edge 

V  Right  Edge  Fracture  Strain 


Fracture  Area  (mm^) 


.055  Y  #1 

.055  Y  #2 

0.7975  | 

Sample 


.035  X  #1 


.035  X  #2 


Left  Edge 

Left  Edge  Fracture  Strain 


1.5995 


1.5019 


.035  Y  #1 


.035  Y  #2 


0.7980 


0.7666 


1.5995 


1.4549 


Fracture  Area  (mm*) 


1.1493 


1.2180 


1.1200 


1.1768 


Table  E3:  Middle  Fracture  Strain  for  Uniaxial  Tests 


# 


73 


Table  E4:  Estimation  of  Average  Uniaxial  Fracture  Stress 


Computed  Values 


■;Av 


Sample 


Edge  Fracture  Strain  (EpsilonE) 


.035  X#1 


.035  X  #2 


1.705 


1.602 


Total  Edge  Area  (AE) 


2.1837 


2.3200 


dSigma/dEpsilon 


300 


300 


.035  Y  #1 


1.625 


2.2120 


300 


.035  Y  #2 


1.478 


2.3262 


300 


.045  X  #1 


1.422 


2.9307 


300 


.045  X  #2 


1.369 


3.0782 


300 


.045  Y#1 


1.440 


2.6788 


300 


.045  Y  #2 


1.458 


2.7353 


300 


.055  X  #1 


1.582 


2.8829 


300 


.055  X  #2 


1.550 


2.9475 


300 


.055  Y  #1 


1.597 


2.8896 


300 


.055  Y  #2 


1.613 


2.8646 


300 


mmm 


Sample 


Fracture  Load  (lb) 


Average  Fracture  stress  (MPA) 


Fracture  Stress  Middle  (MPA) 


.035  X  #1 


697.78 


1019.08 


1120.18 


.035  X  #2 


696.04 


953.24 


1048.34 


Table  E5:  .035”  Uniaxial  Sample  Test  Data 


.035"  Thickness  X  Direction 


file  name:  c2x-1.dat 
description:  0  035"  sheet  metal,  "X"  direction 
'"'  test  date:  11/9/1999 
test  mach.;  tnstron  1331 
specimen  #:  C2X-1 
area:  0 .0176 
thickness:  0.0350 
. width:  0.5035 


file  name:  c2x-2.dat 

description:  0  035"  sheet  metal.  "X"  direction 
test  date  11/9/1999 
test  mach.;  tnstron  1331 
specimen  #:  C2X-2 
area:  0.0176 
thickness:  0  0350 
width:  0.5030 


pretest  notes 

data  acquisition  -  300  msec,  22  bits 
channel  1  =  strain:  10%/v 
channel  2=  load:  562  Ibs/v 

sampling  rate  =  0.5  sec  . _ 

load  rate  =  0.2  in/min 

first  data  point  =  zero  load,  second  data  point  =  zero  strain 


C2X-1 

test  notes 

yeild  stress: 

23,300 

ult  stress: 

46,294 

break  point: 

inside  qaqe  area 

.035  Thickness 


pretest  notes 

data  acquisition  =  300  msec,  22  bits1 
channel  1  =  strain:  10%/v 
channel  2  =  load:  562  IbsAr 
sampling  rate  =  0.5  sec 

load  rate  =  0  2  in/min  !  I 

first  data  point  =  zero  load,  second  data  point  =  zero  strain 


C2X-2 

test  notes 

yeild  stress: 

23,300 

ult.  stress: 

45,903 

break  point: 

inside  gage  area 

Y  Direction 


filename:  c2y-1  dat 

description:  0  035"  sheet  metal.  "Y"  direction 
test  date:  11/9/1999 

test  mach.;  Instron  1331 

specimen#:  C2Y-1 

area:  0  0176 

thickness:  0  0350 
width:  0.5040 


file  name:  c2y-2.dat 

description:  0.035"  sheet  metal.  "Y"  direction 
test  date:  11/9/1999 

test  mach.;  Instron  1331 
specimen  #:'  C2Y-2 

area:  0.0177 

thickness:  0  0350 
width:  0.5045 


pretest  notes  j . 

data  acquisition  =  300 
jchannel  1  =  strain:  10%/ 
channel  2  -  load:  562 


sampling  rate  =0.5  sec 
load  rate  =  0.2  in/min 


msec, 

i/v 


.  22  bits 


first  data  point  = 

zero  load,  second 

data  point  =  zero  strain 

C2Y-1 

test  notes 

1 

yeild  stress: 

25,100 

ult.  stress: 

45,176 

break  point 

inside  qaqe  area 

retest  notes 


data  acquisition  =  300  msec, 
channel  1  =  strain:  10%Ar 
channel  2  =  load:  962  ibs/v 


sampling  rate  =  0.5  sec 
load  rate  =  0.2  in/min 

first  data  point  =  zero  ioad,  second  data  point  =  zero  strain 


C2Y^ 

test  notes 

yeild  stress: 

25.100 

ult  stress: 

45,146 

break  point: 

inside  qaqe  area 
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Engineering  Stress  (psi)  Engineering  Stress  (psi) 


0.035  Sheet  Metal,  MX"  Direction 


0  5  10  15  20  25  30  35  40  45  * 

%  Strain 

Figure  E7:  Engineering  Stress  vs.  Strain  Curve  for  .035”  X  Uniaxial  Test 


0.035”  Sheet  Metal,  "Y"  Direction 


5  10  15  20  25  30  35  40  -  45 

%  Strain 


Figure  E8:  Engineering  Stress  vs.  Strain  Curve  for  .035”  Y  Uniaxial  Test 


Table  E6:  .045”  Uniaxial  Sample  Test  Data 


.045“  Thickness  X  Direction 


file  name:  c3x-1  dat 

description:  0  045'  sheet  metal,  "X’  direction 
test  date:  1  iy9/1 999 

test  mach.;  Instron  1331 
specimen  P:  C3X-1 

area:  0.0227 

thickness:  0  0450 
width:  0.5050 


file  name:  c3x-2.dat 

description:  0  045'  sheet  metal,  “X"  direction 
test  date:  1179/1999 

test  mach.;  Instron  1331 
specimen  #:  C3X-2 

area:  0.0227 

thickness:  0.0450 
width:  05050 


pretest  notes 

data  acquisition  =  300  msec.  22  bits 
channel  1  =  strain:  10%/v 
channel  2  =  load:  562  Ibs/v 
sampling  rate  =  0  5  sec 
load  rate  =  0.2  in/min 

first  data  point  =  zero  load,  second  data  point  =  zero  strain 


pretest  notes 

data  acquisition  =  300  msec,  22  bits  . 
channel  1  =  strain:  10%/v 
channel  2=  load:  562  Ibs/v 
sampling  rate  =  0.5  sec 
load  rate  =  0.2  in/min 

first  data  point  =  zero  load,  second  data  point  =  zero  sir; 


OX-1 

test  notes 

yeild  stress: 

28,600 

ult.  stress: 

47,776 

break  point: 

inside  gage  area 

C3X2 

test  notes 

yeild  stress: 

28,600 

utt  stress: 

47.688 

break  point: 

inside  gage  area 

.045"  Thickness  V  Direction 


file  name:  c3y-1.dat 

description:  0.045"  sheet  metal,  'Y“  direction 
test  date:  11/9/1999 
test  mach.;  Instron  1331 
specimen  #:  C3Y-1 
area:  0.0226 
thickness:  0.0450 
width:  0.5030 


file  name:  c3y-2.dat 

description:  0  045"  sheet  metal.  "Y"  direction 
test  date:  11/9/1999 
test  mach.;  Instron  1331 
specimen  #:  C3Y-2 
area:  0.0226 
thickness:  0.0450 
width:  0.5030 


pretest  notes 

data  acquisition  =  300  msec,  22  bits 
channel  1  =  strain:  10%/v 
channel  2  =  load:  562  Ibs/v 
sampling  rate  =  0  5  sec 
load  rate  =  0.2  in/min 

first  data  point  =  zero  load,  second  data  point  =  zero  strain 


C3Y-1 

test  notes 

yeild  stress: 

30,000 

ult.  stress: 

46,948 

break  point: 

inside  qaqe  area 

pretest  notes 

data  acquisition  -  300  msec.  22  bits; 
channel  1  =  strain:  10%/v 
channel  2  =  load:  562  Ibs/v 
sampling  rate  -  0.5  sec 
load  rate  =  0  2  in/min 

first  data  point  =  zero  load,  second  data  point  =  zero 


C3Y-2 

test  notes 

yeild  stress: 

30,000 

ult.  stress: 

46,999 

break  point: 

inside  gage  area 

77 


Engineering  Stress  (psi) 


0.045”  Sheet  Metal,  "X"  Direction 


Figure  E9:  Engineering  Stress  vs.  Strain  Curve  for  .045”  X  Uniaxial  Test 


0.045"  Sheet  Metal,  "Y"  Direction 


Table  E7:  .055”  Uniaxial  Sample  Test  Data 


file  name:  c4x-1  dal 

description:  0  055“  sheet  metal.  X  direction 
test  date:  6/22/2002 

test  mach.;  UTS 

specimen  #:  C3X-1 

area:  0  0278 

thickness:  0  0550 
width:  0.5050 


.055”  Thickness  X  Direction 

file  name:  c4x-2  dat 

description:  0  055“  sheet  metal.  X  direction 
test  date:  6/22/2002 
test  mach.;  UTS 
specimen  #:  C3X-2 
area:  0  0278 
thickness:  0  0550 
width:  0  5050 


pretest  notes 

data  acquisition  =  300  msec.  22  bits 
channel  1  =  strain:  10%/V 
channel  2  =  load:  562  IbsAr 
sampling  rate  =  0.5  sec 
load  rate  =  0.2  in/min 

first  data  point  -  zero  load,  second  data  point  -  zero  strain 


pretest  notes 

data  acquisition  =  300  msec,  22  bits 
channel  1  =  si  rain:  10%/v 
channel  2  =  load.  562  IbsAr 
sampling  rate  =  0  5  sec 
load  rate  =  0.2  in/min 

first  data  point  =  zero  load,  second  data  point  =  zero  strain 


C3X-1 

test  notes 

yeild  stress: 

8 

oo 

R 

ult.  stress: 

48,202 

break  point: 

inside  gage  area 

OX^ 

test  notes 

yeild  stress: 

28.900 

ult.  stress: 

48,321 

break  point: 

inside  gage  area 

.055”  Thickness  Y  Direction 


file  name:  C4Y-1  dat 

description:  0  055“  sheet  metal,  *Y“  direction 
test  date:  6/22/2002 
test  mach.;  UTS 
specimen  #:  C4Y-1 
area:  0.0277 
thickness:  0  0550 
width:  0  5040 


pretest  notes 

data  acquisition  =  300  msec,  22  bits 
channel  1  =  strain:  10%Ar 
channel  2  =  load:  562  IbsAr 
sampling  rate  =  0.5  sec 
load  rate  =  0.2  in/min 

first  data  point  =  zero  load,  second  data  point  =  zero  strain 


C4Y-1 

test  notes 

yeild  stress: 

30, tOO 

ult  stress: 

47,129 

break  point: 

inside  gage  area 

file  name:  C4Y-2.dat 

description:  0.055’  sheet  metal,  "Y“  direction 
test  date:  6/22/2002 
test  mach.;  UTS 
specimen  U:  C4Y-2 
area: 0  0277 
thickness:  0.0550 
width:  0  5040 

pretest  notes  ,  ] 

data  acquisition  =  300  msec.  22  bits 

channel  1  -  strain:  10%/v 

channel  2  =  load:  562  lbs /v 

sampling  rate  =  0.5  sec 

load  rate  =  0.2  in/min 

first  data  point  =  zero  load,  second  data  point  =  zero  strain 


C4Y.2 

test  notes 

yeild  stress: 

30.215 

ult.  stress: 

47,111 

break  point: 

inside  gage  area 

Appendix  F:  Analytic  Solution  Worksheet 
.035"  Thickness  using  20  mm  Radius  Punch 


Entering  Arguments  for  Analytic  Solution 

Work  Hardening  Exponent  and  Coefficient: 

n  :=  .22  Co:=460106Pa 


Material  Yield  and  Ultimate  Strength 


ay  :=  23300psi  t?u  :=  46100psi 


Plate  Thickness 


tQ  :=  .035in 


Plate  Clamped  Radius  R j  \$mm 


Punch  Radius 


=  20  mm 


Compute  Failure  Parameters 

Vc_fail :=  <-957  +  -399  n)- rad  Vc_fail  =  1  045 

4R  <7 

%ail:=  l*41’n  5^ajj  =  38.781mm 


Compute  Stresses 


Flow  Stress 


+  °u) 


cr0  =  2.392x  108Pa 


Generalized  Membrane  Stress 


Plate  Predictions  at  Failure 


R  :=  110 


Pfai!  :=  2  N0  Rb(sin(vc_fail))2_ 


Pfail =  2-308x  IO  N 


Rcfail :=  Rbsin(^c_fail) 


Ref  :=  — = — 

mm 


r<-:=  Ref..  R 


wfail(rf)  :=  - 


Pfaif'^ 


w 


R  + 


,j(rf)2  -  (sin(^c  fail))4,] 


j(R)2  ~  (sin(vt/ c_fail))4 


2-Tt-N 


mm 


Establish  Load  Range  in  5000  N  Increments,  Up  to  the  Predicted  Failure  Load: 
Load: 

i:=  1..7 


P  := 


0  > 

po 

:=  ON 

Pfail 

5-N 

Pfail 

P1 

5 

pi  = 

=  4.61 6x 

10  N 

Pfail 

4-N 

P2 

Pfail 

4 

p2= 

=  5.771  x 

io3n 

Pfail 

i 

3-N 

•N  P3 

Pfail 

3 

P3  = 

=  7.694  x 

io3n 

Pfail 

2*N 

P4 

Pfail 

P4  = 

=  1.154x 

io4n 

2Pfail 

Z 

3-N 

Pfail 

P5 

_  2Pfail 

3 

P5  = 

=  1.539x 

io4n 

V  N  / 


Wrapping  Angle 

(  0 
22.753 
25.62 

ij;  j  =  |  29.954  deg 
37.699 
44.92 
59.861 ) 


\|/j  :=  asin 


(  I  p  ^ 


2-"-N0-Rb 


\|/0  :=  asin 


vj/2  :=  asin 


\|/4  :=  asin 


2'7t'N0-Rb 


2-n-No-Rb 


2-«-N0-Rb 


\(/6  :=  asin 


fail 
1  2-7t- N0Rb 


v|/l  :=  asin 


\[/3  :=  asin 


2*N0Rb 


2'7t'NC)Rb 


vg 5  :=  asin 


2’tN0.Rb 


Strain  Components 


Err:=  ,n 


ett  '=  -£rr 


zn=  0.143 


-0.103 

-0.143 

-0.234 


-0.689 


Deformed  Plate  Thickness 


t  :=  t0-cosf  vpjj 


True  Material  Stress: 


°rr:=Co(-,n(C0S(vtyi)))n 


2.646x  10 


2.793  x  10 


on=\  3x  10  Pa 


3.342  x  10 


3.64  x  10 


4.238x  10 
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Strain  Components 


zn  :=  In| 


1 


5M 


V  C0SV  '  V  J 


rr 


£tt £rr 


f  0  > 

f  0 

0.081 

-0.081 

0.103 

-0.103 

0.143 

stt  “ 

-0.143 

0.234 

-0.234 

0.345 

-0.345 

^-0.689; 

Deformed  Plate  Thickness 


True  Material  Stress: 


t  :=t0-cos(\j/j) 

arr  “  Co  ‘  (  ln(cos  (  V  i)))n 


arr  ~ 


2.64 6x  10 


,8 


2.793  x  10 


8 


3  x  10 


8 


3.342x  10 


3.64  x  10 


8 


4.238x  10 


.8 


Pa 
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Solve  for  Shape  of  plate,  valid  for  1^.  <  r  <  R 

R:=  110 


0 

4.616 

5.771 

7.694 

11.541 

15.388 

23.082 


N  103 


0 

0.397 

0.447 

0.523 

0.658 

0.784 

1.045 


0 

22.753 

25.62 

29.954 

37.699 

44.92 

59.861 


deg 


RcO:»Rb.Si“(v0) 

mm 

RcO  =  0 

Rcl:=Rb"n('l’l) 

mm 

Rcl  =  7.735 

# 

Rc2»Rt,sin(v2) 

mm 

Rc2  =  8.648 

Rc3:.Rk.iin<'f3) 

mm 

Rc3  =  9.986 

# 

Rc4:=Rh.Sin('l'4) 

mm 

Rc4=  12.23 

Rc5:=Rb.Sin('1,3) 

mm 

Rc5  =  14.122 

• 
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Computation  of  Plate  Shape  at  Several  Loads 


♦ 


Pj  =  4.616x  103N 
r  :=  Rcl..  1 10 


P2  =  5.771  x  103N 
r2:=  Rc2..  110 


P3  =  7.694x  103N 
r3:=Rc3..  110 


4 

P4=  1.154x  10  N 
r4  :=  Rc4„  110 


wl(r):= 


w2(r2)  := 


w3(r3)  := 


w4(r4) 


Pjln 

Lr + 

a 

W  -  (si»(Vl))4l 

-LR  + 

A 

^(R)2  -  (sin(xpl))4_ 

j] 

[  2-n-No 

mm 


— 

ho 

_ El _ 

Lr2+ 

./ J2)2  -  (sin(vp2))4J 

]11 

1r  + 

_J(R)2  -  (sin(\|/2))4] 

I 

2'7t'N0 

mm 


Pjln 

r3  + 

A 

(r3)2  -  (sin(\p3))4_ 

]T1 

R  + 

A 

f(R)2  -  (sin(\j/3))4] 

IJ 

2'7t'N0 

mm 


P4ln 

_r4  + 

J(r4)2  -  (sin 

1 - 1 

]11 

R  + 

A(  R)2-(sin< 

W)"] 

IJ 

2.*-N0 

mm 


P5=  1.539x  104N 
r5:=Rc5..  110 


w5(r5)  := 


P5ln 

r5  + 

U 

(r5)2  -  (sin(v|/5))4]_ 

[r  + 

[J 

^(R)2  -  (sinCvyS))4 _ 

2-n-N0 

mm 


# 
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Displacement  (mm) 


Plot  of  Analytic  Solution  for  Plate  Shape  Up  to  Failure  (Valid  for  r  >  $ 


r>r2,r3,r4,r5,r^ 
Plate  Radius  (mm) 


PUNCH  SIZE: 


Rjj  =  20mm 


r  0  > 

f  0  > 

4.616 

7.735 

5.771 

8.648 

7.694 

N  103  1^  = 

9.986 

11.541 

12.23 

15.388 

14.122 

^23.082> 

v.  17.296> 

Pfail=  2.308  x  104N 


Rc  fail  =  1 7.296mm 


Force  (KN) 


Plot  the  Force  vs.  Displacement  Relation: 


Displacement  at  Point  C: 


Total  Punch  Displacement: 


8: 


°  1 

wl(Rcl) 

w2(Rc2) 

w3(Rc3) 

w4(Rc4) 

w5(Rc5) 


^  wfai](Rc0  J 


5plot :=  % 


(l  -^(Vi)) 


-5; 


mm 


r  0  > 

9.499 

11.477 

14.635 

20.605 

26.311 

^37.637; 


6max:=trunc 


V 


( 1  -  cos  W) 


mm 


-  wfail<Rcf) 


Analytic  Force  vs.  Displacement  Curve 
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045"  Thickness  using  50  mm  Radius  Punch 


Entering  Arguments  for  Analytic  Solution 

Work  Hardening  Exponent  and  Coefficient: 

n  :=  .22  CQ  :=  460106-Pa 

Material  Yield  and  Ultimate  Strength  Gy  :=  28600psi  au  :=  47600psi 


Plate  Thickness 
Plate  Clamped  Radius 
Punch  Radius 


tQ  :=  ,045-in 
R  :=  1 10mm 


Rjj  =  50  mm 


Compute  Failure  Parameters 


Vc  fail :=  (-957  +  -399 n)- rad 


H7 c  fail  ^ 


„  ...  .33  .48  .52 

5fail:=L41n  R  Rb 


Compute  Stresses 


Flow  Stress 


6fail  =  62.453mm 


+  °u) 


a0  =  2.627  x  10  Pa 


^°°7> 


5  N 

N  =  3.467x  10  - 
rn 


Generalized  Membrane  Stress 


Plate  Predictions  at  Failure 


R  :=  )  10 


Pfail:==  .2KNoRb(sin(v('c_fail))2 


Pfai|=  8.l46x  ]0  N 


Rc_fail:=Rbsin('t'c_fail) 


R, 


Ref  := 


c  fail 


mm 


rr  :=  Ref..  R 


& 


w 


fai 


K) := 


Pfail'N 


ff  + 


R  + 


J(rf)2  -  (sin(^c  fail))4 


,j(R)2  ~  (sin(^c_fail))4 


2-rt-N 


o 


mm 


Establish  Load  Range  in  5000  N  Increments,  Up  to  the  Predicted  Failure  Load: 


Load: 

i  :=  1 ..  7 


P  := 


r  n  \ 


Mail 

5-N 

Pfail 

4N 

Pfail 

3N 

Pfail 

2-N 


2P 


fail 


3N 
Pfail 

V  N  J 


N 


P0:=0N 


Pi  := 


rfail 


P2 := 


P3  := 


P4  := 


P5  := 


Pfail 


Pfail 


Pfail 


2Pfail 


P,  =  1.629x  10  N 


P2  =  2.037x  10  N 


P3  =  2.715x  10  N 


P4  =  4.073 x  10  N 


P5  =  5.431  x  10  N 
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Wrapping  Angle 


( 


Vi  = 


f  0 

22.753 
25.62 
29.954 
37.699 
44.92 
V  59.861 ) 


deg 


\j/0  :=  asin 


\|/2  :=  asin 


i|/4  :=  asin 


2.7fN0Rb 


2KN0Rb 


2-^N0-Rb 


\\)6  :=  asin 


[  rfail 

2-«'N0-Rb 


vj/j  :=  asm 


2-"’N0-Rb 


vj/l  :=  asin 


v[/3  :=  asin 


/ 

V> 


2-«-N0-Rb 


xp5  :=  asin 


2-nNoRb 


2'7t'N0-Rb 


m 


Strain  Components 

t  1 


err:=  lnl 


V 


cos^i) 


err  = 


ett  Err 


'  0  > 

(  0  > 

0.081 

-0.081 

0.103 

-0.103 

0.143 

ett  = 

-0.143 

0.234 

-0.234 

0.345 

-0.345 

,0.689; 

v -0.689; 

Deformed  Plate  Thickness 

True  Material  Stress: 


t  :=  t0  cos(v|/j) 
°rr:=Co(~In(cos(vi)))n 


m 


°rr  = 


2.646x  10 


2.793  x  10 


3  x  10 


.8 


3.342  x  10 


3.64  x  10 


.8 


4.238  x  10 


Pa 
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To  Find  the  Maximum  Punch  Displacement: 


R  :=  110mm 


Find  the  Amplitude: 


R 

—  + 

_  Rb 

( R) 
UJ 

-  (sin(v|/j))4 

sin(vj)(l  +  cos^j))  j 

Solve  for  Maximum  Displacement  using  Geometric  Relations 


<:=  -%[l  -  cos(Vi)  +  a-(sin(v|/j))2J 


Find  Radius  of  Contact  Wrapping  Point  C: 


Rc  :=  Rb-sin(xi/i) 


*c  = 


(  0  '' 
19.338 
21.62 
24.965 
30.576 
35.306 
V  43.241 j 


mm 
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Solve  for  Shape  of  plate,  valid  for  <  r  <  R 

R:=  110 


r  0  > 

16.292 

20.365 

P=  27.154  N-103 
40.731 
54.308 
l  81. 462  J 


r  0  > 

22.753 

25.62 

vj/ j  =  29.954  deg 
37.699 
44.92 
i  59.861 ) 


Rc0:=  R.  s'n1'*,0l 
mm 

Rc0=  0 

Rcl^Rb.,i"(*l) 

mm 

Rcl  =  19.338 

Rc2:=Rb.!ln('l'2) 

mm 

Rc2=  21.62 

Rc3sRb.Sin(v3) 

mm 

Rc3  =  24.965 

Rc4:=RbSin(,f4) 

mm 

Rc4  =  30.576 

Rc5  :=  Rl- 

mm 

Rc5  =  35.306 
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Computation  of  Plate  Shape  at  Several  Loads 


P,  =  1 .629 x  ]04N 
r:=  Rcl..  1 10 


wl(r)  := 


P,  In 

Lr+Ll 

(r)2  -  (sin(v|/l))4 _ 

_[r  +  U 

]J 

W  -  (=in(vl))"] 

2^-N0 

mm 


P2  =  2.037x  )04N 
r2  :=  Rc2..  110 


w2(r2)  := 


— 

to 

....  5“ 

_r2  + 

u 

(r2)2  -  (sin(v|/2))4 

]" 

R  + 

J 

11 

^(R)2  -  (sin(v|/2))4_ 

2-k.N0 

mm 


P^  =  2.71 5  x  104N 
r3:=  Rc3..  110 


w3(r3)  := 


_r3  + 

(f3)2  -  (sin(\|/3))4] 

P^-ln 

.  LR  +  N 

f(R)2  -  (sin(\|/3))4]_  _ 

2.n-N0 

mm 


P4  =  4.073  x  104  N 
r4:=Rc4..  110 


w4(r4)  := 


p.-ln 

r4  + 

(r4)2  -  (sin(v|/4))4]_ 

i  ^  in 

.  _R  +  Li 

^(R)2  -  (sln(\^4))4]  _ 

2-k-N0 

mm 


P5  =  5.431  x  104N 
r5  :=Rc5..  110 


P5ln 

[r5  + 

J  (r5)2  -  (sin(\j/5))4]_ 

_Lr  4- 

J( R)2  -  (sin(iy5))4__ 

2.tt.N0 

w5(r5)  := 


Displacement  (mm) 


Plot  of  Analytic  Solution  for  Plate  Shape  Up  to  Failure  (Valid  for  r  >  g): 


Displacement  vs.  Plate  Radius 


PUNCH  SIZE: 


=  50mm 


<  0  > 

r  0  > 

16.292 

19.338 

20.365 

21.62 

27.154 

N-103  1^  = 

24.965 

40.731 

30.576 

54.308 

35.306 

^81.462; 

^43.24L 

Pfail=  8  ,46x  ,()4n 


Rc  fail  =  43.241mm 


Force  (KN) 


Plot  the  Force  vs.  Displacement  Relation: 


Displacement  at  Point  C:  Total  Punch  Displacement: 


f  0  ^ 

f  0  ^ 

wl(Rcl) 

16.893 

w2(Rc2) 

w3(Rc3) 

,  . .  ('-“W)  . 

Splo,  -  Kb'  mm  -«i  6pl«* 

20.125 

5j  ;= 

25.165 

w4(Rc4) 

34.377 

w5(Rc5) 

42.928 

vwfail(Rcf)y 

^59.814; 

8max:=trunc 


%■ 


(l  -  cos(\|/6)) 


-  wfail(Rcf) 


055"  Thickness  using  75  mm  Radius  Punch 


Entering  Arguments  for  Analytic  Solution 

Work  Hardening  Exponent  and  Coefficient: 

n  :=  .22  Co:=460106-Pa 

Material  Yield  and  Ultimate  Strength  CTy  :=  28600 psi  ou  :=47600psi 


Plate  Thickness 
Plate  Clamped  Radius 
Punch  Radius 


tQ  :=  ,055in 
R  :=  1 10  mm 


&  75  mm 


Compute  Failure  Parameters 

Vcjail :=  (-957  +  •399n)  rad  fail  =  1045 

'2'2  40  C9 

5fail  :=  1  41 "  R  V  Sfail  =  77.1 1 1mm 


Compute  Stresses 


Flow  Stress 


) 


XT  •-  O 


a0  =  2.627 x  10  Pa 


Generalized  Membrane  Stress 


Plate  Predictions  at  Failure 


R  :=  110 


Pfail:=_2-7tNoRb(sin(vc_fail))2] 


Pfail=  1  493 x  10  N 


*cjai\:s  Rb*sin(Vc_fei|) 


R, 


Ref  := 


c  fail 


mm 


rf  :=  Ref..  R 


9 


wfail(rf)  := 


_rf + 

[J 

^f)2  "  (sin('t,c_fail))4 

IT 

tail 

R  + 

0 

^R)2  ~  (sin(vP  c_fail))4 

]J 

2-7t-N„ 


mm 


Establish  Load  Range  in  5000  N  Increments,  Up  to  the  Predicted  Failure  Load: 


Load: 

i  :=  1 ..  7 


P  := 


f  c\  \ 


rfail 

5-N 

Pfail 

4-N 

Pfail 

3N 

Pfail 

2- N 

2Pfail 

3- N 

Pfail 

V  ~n"  j 


N 


P0:=0N 


P1  := 


p2:= 


Pfail 


Pfail 


P3  := 


P4 := 


P5  := 


"fail 


"fail 


2Pfail 


Pj  =  2.987x  10  N 


P2  =  3.734x  10  N 


P3  =  4.978x  10  N 


P4  =  7.467x  10  N 


P5  =  9.956x  10  N 
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Wrapping  Angle 


( 


0 


\ 


22.753 
25.62 
29.954 
37.699 
44.92 
V  59.861 ) 


deg 


\[/j  :=  asin 


(  f  P  A 
2nNoRb 


vj/O  :=  asin 


v[/2  :=  asin 


\|/4  :=  asin 


v|/6  :=  asin 


/ 

/ 

U 

f 

V 

/ 

V 


2-*-N0-Rb 


2-TtNoRb 


2-TtNoRb 


|  Mail 
2-™-N0-Rb 


\|/1  :=  asin 


\j/3  :=  asin 


/ 

V 
/ 

V 
/ 


2nN0Rb 


2-w-NoRfc 


\j/5  :=  asin 


2-"-N0-Rb 


Strain  Components 


(  1 

,n  1 

cos 

Vs) 

Ett  ^ 

-£rr 

V  \ 

1  V ) 

'  0  > 

f  0  ' 

0.081 

-0.081 

0.103 

-0.103 

0.143 

8tt  = 

-0.143 

0.234 

-0.234 

0.345 

-0.345 

^  0.689 > 

^-0.689, 

Deformed  Plate  Thickness 

True  Material  Stress: 


t:=  t0cos  (vt/j) 

arr:=Co(-|n(cos(vt/i)))n 


o 


rr 


0  1 

2.646x  10* 

2.793  x  108 

3  x  ]08 

3.342  x  108 

3.64  x  108 

k4.238x  lO8^ 


Pa 


m 


m 


m 
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To  Find  the  Maximum  Punch  Displacement: 


R  1 10mm 


Find  the  Amplitude: 


*|  ef 
_ 1 

-> 

(  R  ^ 

A  2 

-  (sin(vi))4 

1  sin(vj)(l  +  cos(Vi))  I 

Solve  for  Maximum  Displacement  using  Geometric  Relations 


"  cos(vj)  +  a-(sin(\|/j))2] 


Find  Radius  of  Contact  Wrapping  Point  C: 


RC  :=  *b-sin(vi) 


Rc  = 


f  0 

29.007 
32.431 
37.448 
45.864 
52.959 
V  64.861,/ 
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Solve  for  Shape  of  plate,  valid  for  <  r  <  R 

R  :=  110 


P  = 


r  0  > 

'  °  1 

r  0  N 

29.869 

0.397 

22.753 

37.337 

0.447 

25.62 

49.782 

N103  v|/j  = 

0.523 

Vj  = 

29.954 

74.673 

0.658 

37.699 

99.565 

0.784 

44.92 

J  49.347, 

J.045; 

^  59.86 1  j 

deg 


Rc0:=R.-S,nW 

mm 

Rc0=  0 

Kcl»Kb.,in('l',) 

mm 

Rcl  =  29.007 

m 

Rc2:.RhSin(v2) 

mm 

Rc2=  32.431 

Rc3:=Rb5in('l'3) 

mm 

Rc3  =  37.448 

• 

Rc4*RbSin(','4) 

mm 

Rc4  =  45.864 

RcS:=Rb.Sin('l'5) 

mm 

Rc5  =  52.959 

• 

m 


■# 


* 
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Computation  of  Plate  Shape  at  Several  Loads 


mm 


mm 


mm 


mm 
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Displacement  (mm) 


Plot  of  Analytic  Solution  for  Plate  Shape  Up  to  Failure  (Valid  for  r  >  $ 


r,r2,r3,r4,r5,rf 
Plate  Radius  (mm) 


PUNCH  SIZE: 


=  75mm 


f 


0  > 
29.869 
37.337 
49.782 
74.673 
99.565 


N-103 


V  149.347/ 


f  0  ^ 
29.007 
32.431 
37.448 
45.864 
52.959 


mm 


V  64.86lj 


Pfail  =  1  493x  105N 


R,  fajj  =  64.861mm 


Plot  the  Force  vs.  Displacement  Relation: 


Displacement  at  Point  C:  Total  Punch  Displacement: 


f  0  > 

f  o 

wl(Rcl) 

20.79 

w2(Rc2) 

(l  -cos(yj)) 

24.502 

w3(Rc3) 

W-v'  J  -V  «piot  * 

30.165 

w4{Rc4) 

40.193 

w5(Rc5) 

49.228 

vwfai|CRcf), 

v  66.974; 

6max:=trunc 


H- 


(l- 


cos(\|/6)) 

mm 


-  wfail(Rcf) 


Analytic  Force  vs.  Displacement  Curve 


Appendix  G:  Petalling  Work  Approximations 


Petalling  Work  Computation  for  .035”  Sheet 
with  9.5mm  Starter  Hole  and  Four  3-mm  Precracks 


Material  Properties 

Gy  =  23300psi  gu  =  45900psi 


ao  := 


(°y  +  °u) 


Number  of  Petals 


Np:=4 


Crack  Length 


aQ  s  0  mm 


Petal  Semi-Angle 

aj  =  10mm  a5  =  50mm 

=  20mm  a^  =  60mm 

a ^  =  30mm  a y  s  65  mm 

a4  =  40mm  ag  s  70mm 

Moment  Amplification  Factor 

360 

0  := - deg 

2Np 

r|  ^  1 

Plate  Thickness 

Petal  Root  Length 

t  =  .035  in 

b0 :=  ao*sin(0) 

b|  :=  a|-sin(o) 

b^  :=  a^-sin(o) 

b2 :=  a2's*n(0) 

b6  :=  a6sin(0) 

b3  :=  a3  sin(^) 

by  :=  aysin(o) 

b4  :=  a4sin(e) 

b8  :=  a8-sin(e) 

Measured  radius  of  Curvature 


p  :=  150mm 


Bending  Moment 


2V<  b0 

Mo0:=— 7“ 

2a„t2-bi  M  .  2oo-,2b5 
o5=^~“ 

2o|2b,  2c,o  1  b6 

Mo2:=-^  M<*'=— T~ 


2o0-.-b7 


2a  t  bo  A/f  .  "o* 
w  o  j  M„7:= - 

Mo3-= - - -  07  4 


2an  t2b,  M  . 


Starter  Hole  Radius 


rQ  :=  4.75mm 


Cone  Punch  Horizontal  Angle 

<|>  :=  30-deg 


Punch  Vertical  Speed 

V=,a™ 

p  min 

Petal  Perpendicular  Length 


:=  aQ-cos(0) 

:=  ajcos(o) 

l5  :=  a^cos^) 

• 

:=  a2-cos(o) 

16:=  a6cos(e) 

:=  a^-cos^) 

\-j  .=  ay  COS  (9) 

:=  a4cos(e) 

*8  :=a8c°s(e) 

• 

Rate  of  Petal  Length  change 

'dot  :=tan(4.)  v 
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Crack  Tip  Opening  Distance  (CTOD) 
CTO  A 


CTOA  =  10- deg 

5cq  :=  2-(a0-sin(CTOA)) 

8cl :=  2-(arsin(CTOA)) 

8c5  2( 

5c2  :=  2-(a2-sin(CTOA)) 

8c6:=2( 

5c3  :=  2-(aj-sin(CTOA)) 

8c7  :=  2‘( 

8c4  :=  2-(a4-sin(CTOA)j 

8c8  :=  2  ( 

Punch  Travel 

8punch0  :=[tan(e)  (l0+  ro)] 
5punch1  :=[‘an(e)  (l,  +  rj] 
5punch2  :=[tan(e)  (l2  +  rQ)] 
8punch3  :=[tan(e)  (l3+r0)] 
8punch4  :=[tan(0)  (l4+ro)] 


5sin(CTOA)) 

6sin(CTOA)) 

8sin(CTOA)) 

g-sin(CTOA)) 


8punch5  :=[tan(0)  (l5+  ro)] 
8punch6  •=[^tan(0)'('6  +  ro)] 
8punch7  -=  j^tan {()  ( ly  4  rt))j 

5punch8  :=[tan(0)  ('g  +  r0J] 


Rate  of  Work  Done  at  Specified  Increment 


W 


rateO 4 


W 


ratel  = 


Wrate2:= 


W 


rate3 


Mo0, 

"l0-tan(0)> 

( F)  N\ 
°c0 

■~‘dof 

4-n- 

+  3.84 

m 

- 

p  ; 

<  t  j 

Mol  , 

'  lj  tan(0)^ 

f5cl1 

‘dot* 

4rv 

- . 

+  3.84| 

— 

m 

r 

\  9  J 

1  t  ) 

Mo2, 

( l2’tan(0)  N 

fy 

- ’dot' 

4n- 

+  3.84 

— 

m 

r 

V  P  ) 

1  t  J 

Mo3, 

r  lytan(0)^ 

f  y 

m  dot* 

4-n- 

+  3.84| 

— 

m 

_ 

l  p  J 

1  t  J 

-4 


(sin(e))  3  -(cos(e))  1 


—  (sin(e))  -(cos(e)) 


-  ! 


-  4 
3 


-  1 


M 


W 


o4 


rate4  •  'dot 

m 


M 


W 


rate5  * 


o5 

m 


4*n* 


4’t]  * 


+  3.84 


ry 

t  j 


M 


W 


06 


rate6  ■ 


m 


IVL 


**dot 


4n 


W,..e7  =  -f!'do.' 
m 


4-n 


l4-tan(e)  ^ 


P  J 


lrtan(o) 


(  l6»an(0)  ^ 


v  p  ; 


Ijtan^)^  f S^V  /  n3 


(sin(e))  -(cos(e)) 


-4 

f  |  (sin(e))  3  (cos(e))  1 


—  -  4 

'^sin^e^  3  (cos(0)) 


(sin(©))  3  (cos(e)) 


+  3.84 


M 


W 


08 


rate8  •“  ‘dot 


m 


4-tj  ■ 


V  p  J 


lolanCe) 


P  J 


+  3.84 


V  t  J 


(sin(e))  3  (cos(e))  1 


V  t  ) 


-4 

3  /  /^w-  1 


“  I  (sin(e))  *(cos(e)) 
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Total  Membrane  Work 


W 


mO' 


M 


(’o  -  ro) 


oO’ 


m 


3.84 


1 


[t]  (f]3-(sin(9))  3 


Wml:= 


Wm2:= 


wm3:= 


Mo  r 


Ol-o) 


m 


1 


—  t 


3.84 


P 


M 


fe-o 


o2* 


m 


3.84 


M 


o3‘~ 


m 


3.84 


V  »  J  \X  ) 

n3 

t  j 

s  \ 

°c3 

*  J 


(sin(e)) 


^Sin^9^ 


—  L 

3 


1 


•If J  (-(e)) 


w 


m4* 


Wm5:= 


M 


('4  “  ro) 


o4' 


3.84 


Jc4 

t  ; 


—  ]  (sin(e)) 
t 


M 


o5‘~ 


Os-g 


m  , 


1 


3.84 


n3 

^  t  ) 


(f)  <sin(e)) 


wm6~ 


M 


(!6-ro) 


06' 


m 


3.84 


W 


m7‘ 


M0  r 


(*7  “  ro) 


r*_c6\3  '  "3 

U  J 

1 


—  t 


f]  (sin(9» 


3.84 


^)3(£L)3.W9»  3 


Wm8:= 


M 


08* 


Os-g 


m 


3.84 


/  x3 

V  t  j 


•  (sin(G)) 

t 


—  i 

3 


(cos(e))  1 

(cos(e))  1 

(cos(e))  1 

(cos(e))  1 

(cos(e))  1 

(cos(e))  1 

(cos(e))  1 

(cos(e))  1 

(cos(e))  1 
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Total  Bending  Work 


'Vb0-'VMo0('0-ro) 


2  tan(0) 


wbl  -  NP  Mol  ('l  -  ro)2 ~  wbS  »  Np  Mo;.(l5-  r0)2 

wb2  =■  vMb  -  02  ^  Wbs  *  np.Mo6.|,6-  tfJS®. 

Wb3  *  Np  Mo3  ('3  -  'of  ~  »b7  »  Np  Mor('7  ‘  'of  ~ ^ 


wb4*Np'Mo4  0-|-'o): 


2  pm 
2  tan(e) 
2  pm 


wb8-NpMo8(l8-r0) 


2  pm 
2  tan(Q) 
2  pm 


Membrane  and  Bending  Work  Totals 


Wmo=ONm 

Wb0  =  0N-m 

Wm]  =  0.641Nm 

WbJ  =  4.788x  10" 

3Nm 

Wm2=  6.54Nm 

Wb2  =  1.568x  10" 

3Nm 

Wm3=  13.648N  m 

Wb3  =  7.227x  l(f 

3Nm 

Wm4=  41.295Nm 

Wb4  =  0.02N-m 

Wm5=  72.31  INm 

Wb5  =  0.042N-m 

Wm6=  113.5 14Nm 

Wb6  =  0.076Nm 

Wm?=  141.606N  m 

Wb?  =  0.098Nm 

Wm8=  165.581N  m 

Wbg  =  0.125Nm 

Total  Work 

WtO:=Np(WbO+Wmo) 

Wtl:=Np-(Wbl  +  Wml) 

Wt2:=Np(Wb2+Wm^ 

Wt3:=Np-(Wb3+Wm3) 

Wt4:=MWb4+Wm4) 
Wt5:=MWb5  +  Wm5) 

Wt6^=Np(Wb6+Wm6) 

Wt7:=Np(Wb7+Wm7) 

Wt8;=Np(Wb8+Wm8) 


wt0  =  oj 

Wtl  =  2.566J 
Wt2  =  26.167J 
Wt3  =  54.622J 
Wt4=  165.257J 
Wt5  =  289.409J 
Wt6  =  454.36J 
Wt?  =  566.816J 
Wt8  =  662.824J 
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Estimation  of  Forces 


W 


F0:=  2", 


F,  :=  2- 


tO 


^punchO 


W 


tl 


W 


F5:=2*> 


t5 


^unchl  &punch5 


W 


F2:=2“ 


12 


W 


°punch2 
W* 


f6^2: 


t6 


3punch6 


F3  :=  2*" 


t3 


W 


°punch3 


F7 :=  2*“ 


t7 


°punch7 


W 


F4:=2- 


t4 


W 


F8 :=  2" 


t8 


C  0  C 

°punch4  °punch8 


Force  and  Displacemnt  Arrays 


V 

f  0  ^ 

f  0  ^ 

R 

434.109 

/  ~  \ 

F1 

^punchl 

u 

F2 

^punch2 

2.77  x  103 

11.821 

_ _  3 

18.892 

F3 

^punch3 

4.208x  10 

A 

25.963 

F4 

6  := 

c 

°punch4 

P  = 

1.001  X  10 

N  6  = 

33.034 

F5 

c 

°punch5 

1.443x  104 

40.105 

r 

4 

47.176 

Upunch6 

1.926x  10 

50.712 

F7 

^punch7 

2.235  x  104 

! 

^54.247; 

^punch8  ) 

^  2.444  x  104y 

b 
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% 


Appendix  H:  LS-DYNA  Input  File 

$  .055"  plate  with  75  mm  radius  sphere 
$ 

^KEYWORD 

’TITLE 

Stamping 

TART 

part 

1  2  2 

$... + — 1 — - + — 2- — + — 3 — + — 4  — - + — 5 —  +  — 6 — + — 7 — + — 8 
$  (2)  CONTROL  CARDS. 

$ — + — i. — +  -—2- — + — 3 — + — 4 — + — 5 — + — ^ — + — 7 — + — 8 
^define  curve 
$  lcid  (x-direction) 

1 

$  abscissa  ordinate 
0.0  1000.0 
100.0  1000.0 
*control_termination 

$  ENDTIM  ENDCYC  DTMIN  ENDNEG  ENDMAS 
l.e-1 

*control_timestep 

$  DTINIT  SCFT  ISDO  TSLIMT  DTMS  LCTM  ERODE  MS1ST 
.000  .900  0 

*CONTROL_SHELL 

$  WRPANG  ITRIST  IRNXX  ISTUPD  THEORY  BWC  MITER 
20.000  2  -1  1  2  2  1 

*CONTROL_HOURGLASS 
$  IHQ  QH 

1  .100 

’tCONTROL_BULK_VISCOSITY 
$  Q2  Q1 
1.500  .060 

*CONTROL_ENERGY 
$  HGEN  RWEN  SLNTEN  RYLEN 

2  2  11 

TONTROLOUTPUT 

$  NPOPT  NEECHO  NREFUP  IACCOP  OPIFS  IPNINT  IKEDIT 
0  0  0  0  .000  0  100 

TATABASEEXTENTBINARY 

0  0  3  0  1  1  1  1 

0  0  0  0  0  0 

i — 2 — + --3 — + --4  -- + — 5 — + — 6 — + —7— + —8 
$  (4)  DATABASE  CONTROL  CARDS  FOR  ASCII  FILE 

$ — |- — 1 - 1-  — 2 - 1  — 3 — + — 4 — + — 5 — + — 6 — + — 7 — — 8 


% 

% 


* 
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"DAT ABASERWFORC 
l.e-4 

*DATABASE_GLSTAT 

l.e-4 

*DATABASE_MATSUM 

l.e-4 

$ — i- — i f  — 2 1 — 3 1 — 4 v  — 5 — -i — 6 — (-  — 7 — i — 8 

$  (5)  DATABASE  CONTROL  CARDS  FOR  BINARY  FILE 

$-- + -- 1  --+ - — 2 — + —3 — + —4  -- + —5 — + --6— + —7— + —8 
*DATABASE_BINARY_D3PLOT 
$  DT/CYCL  LCDT  NOBEAM 
l.e-3 

$-+ — 1 — + — 2- — + — 3  -- + —4 — + —5 — + -6-- + —7— + — 8 
$  RIGID_SPHERE 

— i — + — 2-—+ — 3 —  +  — 4 — + — 5 — + — 6— +  — 7 —  +  — 8 
*rigidwall_geornetric_sphere motion 
$  NED/NSID 
0 

S  xt  yt  zt  xh  yh  zh  fric 

0.  100.  75.7  0.  100.  0.  0.05 

$  radius 
75.0 

$ 

$  LCID  OPT  VX  VY  VZ 
1  0  0.  0.  -1.0 

$ — i- — i f- — 2 — + — 3 — + — 4 h — 5 — + — 6 — h 7 —  +  — 8 

$  BOUNDARY  CONDITION 


$ - 1 - 1 - 1- 

— 2 — + — 3 - h- 

-4 - | - 5__ 

--+ — 6 — 

-+ — 7 — +■ 

— 8 

^BOUNDARY  SPC  SET 

Sclamped 

$  NID/NSID 

CED  DOFX 

DOFY 

DOFZ 

DOFRX 

DOFRY 

DOFRZ 

1 

1  1  1 

1  1 

1 

^BOUNDARY  SPC  SET 

$symy 

$  NID/NSID 

CED  DOFX 

DOFY 

DOFZ 

DOFRX 

DOFRY 

DOFRZ 

2 

1 

1 

1 

^'BOUNDARY  SPC  SET 

$symx 

$  NID/NSID 

CID  DOFX 

DOFY 

DOFZ 

DOFRX 

DOFRY 

DOFRZ 

3 

1 

1 

1 

"BOUNDARY  SPC  SET 

Scenter 

$  NID/NSID 

CID  DOFX 

DOFY 

DOFZ 

DOFRX 

DOFRY 

DOFRZ 

4 

1  1 

1  1 

1 

$ — i - 1 - 1-. 

— 2 — + — 3— — 

--4 — + — 5- 

-+ — 6 — 

•+ — 7 —  +  - 

-8 

=:MAT_PIECEWISE_LINEAR_PLASTICITY 
$  MID  RO  E  PR  SIGY  ETAN 
2  7.800E-09  210.E+03  0.33 


EPPF  TDEL 


S  C  P  LCSS  LCSR 

$  EPS1  EPS2  EPS3  EPS4  EPS5  EPS6  EPS7  EPS8 

0  .02  .04  .06  .1  .5  .75  1 

$  ESI  ES2  ES3  ES4  ES5  ES6  ES7  ES8 

240  260  295  330  695  920  1160  2000 

*SECTION_SHELL 
SHMNAME  PROPS  lSectShlll 

$  SID  ELFORM  SHRF  NIP  PROPT  QR/IRID  ICOMP  SETYP 
2  2  0.833  5 

$  T1  T2  T3  T4  NLOC 

1.397  1.397  1.397  1.397 

$— + — i — + — 2- — + — 3 — + — 4 — + — 5 — + — 6 — + — 7 — + — 8 
S  GEOMETRY 

$ — h — 1 1 — 2 — i — 3  — i 4  — i — 5  — i — 6 — i 7 i 8 

*NODE 
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Appendix  I:  Strain  Field  Worksheet 


Strain  Field  Computation  Worksheet 

ROW  0 


ui  i_i_ptoo  =  0  un_2_ptoo  =  0 

^1 1  _pt00  -  1  0  §y  22_pt00  =  1  0 

u22_l_pt00  =  0  u22_2_pt00  =  0 

^MljrtOO :=  U1  l_2_pt00  "  U1  l_l_pt00 

8u22_pt00:-  u22_2_pt00 

^lljtOO 

5u22_pt00 

Ell_pt00~  , 

“1  l_pt00 

£22_pt00  -  x 

5y22_pt00 

lIT^H-PtOo)  f  5u22_pt00  Yl 

'  1 2_pt00  • 


2Ll  5y22jpt00  J 


V  5x]  1  _ptOO  JJ 


=-ll_ptOO 

E22_pt00 
E 1 2_pt00 


=  0 
=  0 
=  0 


Next  Point 


• 

ull_l_ptl0=  8 
^1  l_ptl0  =  -8 
u22_l_pt10=  0 

ull_2_ptl0  =  10 
^22_ptl0=  8 
u22_2_ptl0-  0 

• 

Sull_ptl0:-ull_2_ptl0~  ull_l_pt10 

®u22__ptl0 u22_2_j)tl0“  u22_l_ptl0 

^lljtlO 
en_ptio  -  , 

6xll_ptlO 

6u22_ptlO 
£22  ptlO  —  , 

S>22ptl0 

• 

1 

e  1 2_ptl  0 "  2 

<5ul l_ptlO^|  i  ^u22_ptlo'' 
v5y22_ptloJ  V6x]l_ptl0; 

el  l_ptl0  =  0  25 
E22_ptl0  =  0 
E]2_ptl0  =  0125 

123 


Epsilon-n  Strain  Field: 


£ 


11:= 


E1  l_pt06  E1  l_ptl6  e  1 1 _ pt26  £1  l_pt36 

ell_pt05  el  l_ptl5  e  1  l_pt25  E1  l_pt35 
el l_pt04  el l_ptl4  el l_pt24  el l_pt34 
E1  l_pt03  el  l_pt13  E1  l_pt23  E1  l_pt33 
el  l_pt02  e  1 1  __pt1 2  el  l_pt22  e  1 1  _pt32 
Ell_pt01  Ell_ptll  E1  l_pt21  £1  l_pt31 
<Ell_ptOO  Ell_ptlO  el  l_pt20  E1  l_pt30 


E 1  l_pt46  e  1  l_pt56  E1  l_pt66  E1 1  jjt76^ 
el  l_pt45  ell_pt55  Elljpt65  el l_pt75 
el  l_pt44  E1  l_pt54  Ell_pt64  ell_pt74 
E1  l_pt43  E1  l_pt53  ell_pt63  E1  l_pt73 

E1  l_pt42  E 1  l__pt52  Ell_pt62  el  1 _ pt72 

£1  l_pt41  Ell_pt51  el  l_pt61  E1 1 _ pt71 

Ellj3t40  ell_pt50  ell_pt60  E1  l_pt70 > 


f  0 

0.1 

0.15 

0 

0 

0.2 

0.4 

0.25  ^ 

0.5 

0.1 

0.2 

0.05 

0.2 

0.3 

0.4 

0.3 

-0.222  0.111 

0.2 

0.1 

0.2 

0.4 

0.4 

0.5 

E11  - 

-0.25 

0.222 

0.222  0.278  0.333 

0.444 

0.722  0.579 

0 

0.111 

0.167  0.222  0.333 

0.444 

1.056  0.842 

0 

0.111 

0.25 

0.25 

0.5 

0.5 

0.625 

0.941 

<  o 

0.25 

0.25 

0.375  0.571 

0.75 

0.5 

0.75  , 

Epsilon??  Strain  Field: 


£22 :" 


s22_pt06 

£22_pt05 

£22_pt04 

s22_pt03 

s22_pt02 

£22_pt01 

V£22_pt00 


£22_ptl6 
e22_ptl  5 
£22_ptl4 
e22_ptl3 
£22_ptl2 

£22_ptl  1 
£ 


e22_pt26 

e22_pt25 

E22_pt24 

e22_pt23 

£22_pt22 

£22_pt2 1 
E' 


22_ptl0  t22_pt20 


£22_pt36 
£22_pt35 
E22_pt34 
£22_pt33 
e22_pt32 
£22_pt3 1 
£22_pt30 


e22_pt46 

£22_pt45 

£22_pt44 

e22_pt43 

e22_pt42 

E22_pt41 


t22_pt40  t22_pt50 


E22_pt56 

e22_pt55 

e22_pt54 

e22_pt53 

£22_pt52 

e22_pt5 1 
£' 


e22_pt66 

E22_pt65 

e22_pt64 

e22_pt63 

E22_pt62 

E22_pt61 

E22_pt60 


e22_pt76>1 

e22_pt75 

£22_pt74 

e22_pt73 

£22_pt72 

e22_pt71 

£22_pt70  J 


z22~ 


0.4 

0.5 

0.55 

0.6 

0.7 

0.7 

0.8 

0.8 

0 

0.5 

0.55 

0.65 

0.7 

0.9 

1 

3 

0.2 

0.3 

0.45 

0.45 

1.7 

1 

1 

1.05 

0.2 

0.111 

0.222 

0.333 

0.389 

0.944 

1 

1 

0 

0 

0.167 

0.222 

0.222 

0.222 

3.1 

3.15 

0 

0 

0 

0 

0.1 

0.15 

0.2 

0.25 

0 

0 

0 

0 

0 

0.125 

0.187 

0.25 

125 


Epsilon1?  Strain  Field: 


£12:~ 


<e12_pt06  £12_ptl6  £  1 2_pt26  £12_pt36  e12_pt46  e12__pt56 
£  1 2_pt05  £12_pt15  £12_pt25  £12_pt35  £12_pt45  e12_pt55 
c12_pt04  £12_ptl4  e  1 2_pt24  e12_pt34  e  1 2_pt44  e12_pt54 
£  1 2_pt03  E 1 2  pt  1 3  e  1 2_pt23  £12_pt33  £12_pt43  e12_pt53 
£  1 2_pt02  e12_ptl2  e  1 2_pt22  £12_pt32  £12_pt42  £12_pt52 
£12_pt01  £12_ptll  e12_pt21  e  1 2_pt3 1  e12_pt4]  e12_pt51 
^ e  1 2_pt00  e  1 2_pt  1 0  e12_pt20  e12_pt30  £12jpt40  £12_pt50 


£12_pt66 

£12_pt65 
E12_pt64 
e  1 2_pt63 

e12_pt62 

e12_pt61 

e12_pt60 


f  0.2 

0.3 

0.35 

0.3 

0.35 

0.45 

0.6 

0.525^ 

0.25 

0.3 

0.375 

0.35 

0.45 

0.6 

0.7 

0.65 

0.011 

0.217 

0.325 

0.275 

0.95 

0.7 

0.7 

0.775 

0.025 

0.167 

0.222 

0.306 

0.361 

0.694 

0.881 

0.801 

0 

0.063 

0.167 

0.222 

0.278 

0.333 

1.086 

1.005 

0 

0.063 

0.125 

0.125 

0.263 

0.294 

0.375 

0.547 

<  o 

0.125 

0.125 

0.187 

0.25 

0.437 

0.344 

0.5  > 

Rotating  the  strain  field  in  order  to  align  coordinate  sysytem  with  crack: 

e22 

£rm  £22'cos<22  de8)  £^:=  — (22  deg)  £4n  :=  £12’C°s(22  deg) 


£12_pt76> 

E 1 2_pt75 

£  1 2_pt74 

e12_pt73 

e12_pt72 

e12_pt71 

e12_pt70> 


Final  Crack-Oriented  Strain  Fields: 


-Tin 


'  0.371  0.464  0.51  0.556  0.649  0.649  0.742  0.742^ 
0  0.464  0.51  0.603  0.649  0.834  0.927  0.927 

0.185  0.278  0.417  0.417  1.576  0.927  0.927  0.974 
0.185  0.103  0.206  0.309  0.361  0.876  0.927  0.927 
0  0  0.155  0.206  0.206  0.206  1.02  1.066 

0  0  0  0  0.093  0.139  0.185  0.232 

<  0  0  0  0  0  0.116  0.174  0.232 \j 


Ett  = 


( 0.431  0.539  0.593  0.647  0.755  0.755  0.863  0.863^ 

0  0.539  0.593  0.701  0.755  0.971  1.079  1.079 

0.216  0.324  0.485  0.485  1.834  1.079  1.079  1.132 

0.216  0.12  0.24  0.36  0.419  1.019  1.079  1.079 

0  0  0.18  0.24  0.24  0.24  1.186  1.24 

0  0  0  0  0.108  0.162  0.216  0.27 

K.  0  0  0  0  0  0.135  0.202  0.27 , 


e4n  = 


f  0.1 85  0.278  0.325  0.278  0.325  0.417  0.556  0.487^ 

0.232  0.278  0.348  0.325  0.417  0.556  0.649  0.603 

0.01  0.201  0.301  0.255  0.881  0.649  0.649  0.719 

0.023  0.155  0.206  0.283  0.335  0.644  0.816  0.743 

0  0.058  0.155  0.206  0.258  0.309  1.007  0.932 

0  0.058  0.116  0.116  0.243  0.272  0.348  0.507 

<  0  0.116  0.116  0.174  0.232  0.406  0.319  0.464, 
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